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00 ; ABSTRACT 

p\| , We present the initial imaging and spectroscopic data acquired as part of the VLT VIMOS 

^^ ■ Lyman-break galaxy Survey. UBR (or UBVI) imaging covers five w 36' x 36' fields centred 

CO ' on bright z > 3 QSOs, allowing w 21, 000 2 < z < 3.5 galaxy candidates to be selected 

(/-^ . using the Lyman-break technique. We performed spectroscopic follow-up using VLT VIMOS, 

(•~^ ■ measuring redshifts for 1020 z > 2 Lyman-break galaxies and 10 z > 2 QSOs from a total 

^^ I of 19 VIMOS pointings. From the galaxy spectra, we observe a 625 ± 510 kms^ velocity 

offset between the interstellar absorption and Lya emission line redshifts, consistent with 
previous results. Using the photometric and spectroscopic catalogues, we have analysed the 
galaxy clustering at z « 3. The angular correlation function, w{9), is well fit by a double 
^\f ' power-law with clustering scale-length, tq — 3.19^0 54 /i^^Mpc and slope 7 = 2.45 for 

^ . r < 1 h^^Mpc and tq = 4. 37^0:55 h'^Mpc with 7 = 1.61 ± 0.15 at larger scales. Using the 

redshift sample we estimate the semi-projected correlation function, Wp{a) and, for a 7 = 1.8 
power-law, find r^ = ?,.Qltlfi /i^^Mpc for the VLT sample and tq = i-9^tl\t /i-^Mpc 
for a combined VLT+Keck sample. From ^(s) and ^(u, tt), and assuming the above £,(r) 
models, we find that the combined VLT and Keck surveys require a galaxy pairwise velocity 
dispersion of pa 700 kms^ , higher than the « 400 kms^ assumed by previous authors. We 
also measure a value for the gravitational growth rate parameter of /3(z = 3) = 0.48 ± 0.17, 
again higher than previously found and implying a low value for the bias of fe = 2.061q 5. 
This value is consistent with the galaxy clustering amplitude which gives b = 2.22 ± 0.16, 
assuming the standard cosmology, implying that the evolution of the gravitational growth rate 
is also consistent with Einstein gravity. Finally, we have compared our Lyman-break galaxy 
clustering amplitudes with lower redshift measurements and find that the clustering strength 
is not inconsistent with that of low-redshift L* spirals for simple 'long-lived' galaxy models. 

Key words: galaxies: intergalactic medium - kinematics and dynamics - cosmology: obser- 
vations - large-scale structure of Universe 
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1 INTRODUCTION 

Observations of the 2 ~ 3 galaxy population present a valuable tool 
for studying cosmology and galaxy formation and evolution. For 
cosmology, the interest is in measuring the galaxy clustering am- 
plitudes and redshift space distortions at high redshift. They both 
lead to virtually independent estimates of the bias whose consis- 
tency leads to a test of the standard cosmological model. For the- 
ories of galaxy formation and evolution, this is a key period in the 
history of the Universe in which significant levels of star formation 
shape both galaxies and the inter-galactic medium (IGM) around 
them. An especially vital direction of study is the effect of galac- 
tic winds at this epoch. Such winds have been directly obse rved a t 
low JHeckman et al.lll990l:lLe"hnert et alJl999l.lM artin 20 ol.l2006l) 
and high ('Pettini et al."200l';'Adelb erger et al]|2003l : IWilman et al.1 
12005 ; Adelberger et al. 2005a) redshift and are invoked to explain 
a range of astrophysical phenomena. 

A basic item of cosmological interest is the spatial cluster- 
ing of the z ~ 3 galaxy population itself. In ACDM, structure 
in the Universe is known to grow hi erarch i cally through grav- 
itational instability (e.g. .Mo & Whit3 1 19961 ; Ijenkins et al] 1 19981 : 
ISpringel et al.ll2006n and testing this model requires the measure- 
ment of t he clusteri ng of mat ter in the Universe across cosmic 
time (e.g. ISpringeTe t al. 2005; IOrsi et"ak 2008; Kim et al. 2009). 
Surveys of matter at z ~ 3 currently focus on two main popu- 
lations, LBGs and Lyman-a emitters (LAEs). A number of mea- 
sureme nts of galaxy clu s tering are available at z ~ 3. Fo r ex- 
ample, lAdelberger et al.l J2003f) and Ida Angela et al. (2005a) use 
the Ke ck LBG sample with spectroscopic redshfits of Steidel e t al.l 
( 120031) to measure LBG clustering clustering lengths of ro = 
3.96 ± 0.29 /i-^Mpc and ro = 4.48l^:i^ /i"^Mpc respectively. 
Further surveys of LBGs at 2: ~ 3 have pr oduced a range of results 
with, for example. iFoucaud et al. (2003) measuring a clustering 
length for a photometric sample s elected from the CFHT L egacy 
Survey of ro = 5.9±0.5 /i~^Mpc. lAdelberger et Zl ( l2005bl) mea- 
sured ro = 4.0 ±0.6 /i~^Mpcat (z) = 2.9 us ing a different photo- 
metric sample whilst lHildebrandt et al.l <2007h measured a value of 
ro = 4.8±0.3 /i^^Mpc from an LBG sample taken from GaBoDS 
data. 

Ida Angela et al.l j2005an go on to use the Keck LBG sam- 
ple to investigate, via redshift space distortions, the gravitational 
growth rate of the galaxy population at 2: ~ 3, measuring an in- 
fall parameter of /3(z = 3) = 0.25^006- The infall parameter, 
/3, quantifies the l arge-scale infall towa rds density inhomogeneities 
( lHamilton|[T992l ; iHawkins et al.ll2003l) and is defined as I3(z) = 
0,m{z)^'^ /h{z), where Q,m{z) is the matter density and 6(2) is 
the bias of the galaxy population. The 2dF Galaxy Redshift Survey 
(2dFGRS) measurement of the infall param eter nearer the presen t 
epoch gave P{z « 0.1) = 0.49 ± 0.09 dHawkins et ai]|2003h . 
similar to values ob tained by previous local measurements (e.g. 
iRatcliffe et alj|l998r) . There have also been dynamical measure- 
ments of P at interm ediate redshifts using Luminous Red Galax- 
ies where IRoss et al. (2007) found P{z = 0.55) = 0.4 ± 0.05. 
Ida Angela etal., C2005b) used the combined 2dF and 2SLAQ QSO 
redshift surveys to find P{z = 1.5) = 0.60 ± 0.14. Finally, 
IGuzzo et al.l ( 120080 used the VVDS galaxy redshift su rvey to mea- 
sure P{z = 0.77) = 0.70 ± 0.26. As emphasised bv lGuzzo etall 
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( l2008f) . if there are independent estimates of b{z) for each red- 
shift sample, then the standard model prediction for the evolution 
with redshift of the gravitational growth rate of / — Q,m (2)"" can 
be tested agains t alterna tiv e gravity models. H ere, we shall follow 
Ida Angela et all ( l2005allbl) ; lHovle et alj ( l2002h in making their ver- 
sion of the reds hift- space distortion cosmol ogical test which also 
incorporates the JAlcock & Paczvnskil j 19791) geometric cosmolog- 
ical test. 

From redshift-space distortions, we can also determine the 
small-scale dynamics of the galaxy population which are usually 
simply modelled as a Gaussian velocity dispersion, measured from 
the length of the 'fingers-of-God' (Jackson 1972; Kaiser 1987|) 
in redshift space clustering. This velocity dispersion will gener- 
ally als o include the effects of velocity measurement error. Al- 
though Ida Angela et al.l ( l2005ah had to assume a fixed value of 
< wl >i'2= 400 kms~^ for the mean pairwise velocity disper- 
sion when making their LBG measurement of /?(2 = 3), in bigger 
surveys it is possble to fit for < wl, >^" an d /3 sim ultaneously. 
Thus in 2dFGRS at 2 « 0.1. iHawkins et al.l ( l2003h measured a 
pairwise velocity dispersion of < wl >^'^« 500 kms~^. As well 
as being of interest cosmologically, the intrinsic galaxy-galaxy ve- 
locity dispersion is interesting in terms of establishing the group en- 
vironment for galaxy formation. Furthermore, these random pecu- 
liar velocities dominate at the smallest spatial scales, significantly 
affecting clustering measurements on scales r < 5 /i~^Mpc. They 
influence both the observed galaxy-galaxy clustering and the ob- 
served correlation between galaxy position s and nearby Lyg forest 
absorp tion from the IGM (as m easured in lAdelberger et alj|2003l 
l2005al and lCrighton et al.ll20ldl) . To interpret galaxy-IGM cluster- 
ing results we shall see that measurements of the small scale dy- 
namical velocity dispersion of the galaxy population are very im- 
portant. 

Galactic winds powered b y supernovae are a crucial ingredient 
in models of galaxy formation (^ Dekel & Silkll98q ; lwhite & FrenkI 
L199I). Such negative feedback is required to quench the forma- 
tion of small galaxies and make the observed faint-end of the 
galaxy luminosity function much flatter than the low-mass end of 
the dark- matter mas s funct ion, see for example the semi-analytical 
model of ICoIe et alj J2OO0l) . Simulations without such strong feed- 
back tend to produce galaxies with too massive a bulge, which 
consequently do not lie on the observed Tully-F isher relation 
dSteinmetz & Navarro! 1 19991 ; iGovemato et al.l I2OI0I) . Such winds 
can also remove a significant fraction of baryons from the forming 
galaxy, t hereby explain ing why galaxies are missing most of their 
baryons jBregman et al. 2009). and hence a re much fainter in X-ray 
emission than expected dCrain et al.ll2010l) . In addition, observa- 
tions of the IGM as probed with QSO sightlines reveal the presence 
of metals even in th e low density regions producing Lyg forest ab- 
sorption dSongaila & Cowiell996l ; |Pettim et al.l2003l ; lAguirre et all 
l2004l; lAraciI et alj2004h . Other than enrichment from galactic scale 
winds, it is difficult to see from whe re these metals origi nate and 
this is confirmed by simulations (e.g. lWiersma et alj 20091) . 

Direct evidence for outflows in high redsh ift galaxies came 
from t he K eck LBG surv e y spec tra analysed bv lAdelberger et al.l 
d2003r) and Ishaplev et al.l d2003r) who found evidence for offsets 
in the positions of ISM absorption lin es, Lyg emission and rest - 
frame optic al emission lines (see also IPettini et alJl200CI . 120020 . 
IShaplev et a l. (2003) present a model in which the optical emis- 
sion lines arise in nebular star-forming HII regions, giving the in- 
trinsic galaxy redshift, whilst the ISM absorption lines originate 
from outflowing material surrounding the stellar/nebular compo- 
nent. Lyg emission arises in the stellar component, but outflowing 
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neutral material scatters and absorbs the blue Lya wing, leaving 
a peak redshifted with respect to the intrisic galaxy redshift (e.g. 
ISteidel et alJl20 10). One of our prime aims here is to test the ob- 
servations underpinning this model in an independent sample of 
LBGs. 

In this paper, we present the first instalment of data of a 2 --^ 3 
survey of LBGs within wide (~ 30') fields centred on bright z ^ ?, 
QSOs. We discuss the imaging and spectroscopic observations, the 
latter including a search for redshift offsets in the LBG spectra, 
followed by an analysis of the clustering and dynam ics of the LBG 
galax y populations in our fields. In a further paper JCrighton et al.l 
|20IQ !). we present the analysis of the relationship between LBGs 
and the surrounding IGM via QSO sight-lines, with the intent of 
further investigating the extent and impact of galactic winds on the 
IGM. 

The structure of this paper is as follows. We provide the details 
of our imaging survey in section |2] covering observations and data 
reduction. In section |3] we present VLT VIMOS spectroscopic ob- 
servations, describing the data reduction and object identification 
processes. Section|4]presents a clustering analysis of the photomet- 
rically and spectroscopically identified objects and we finish with 
our conclusions and summary in section[5] Unless stated otherwise, 
we use an n^ = 0.3, Q,a = 0.7, Ho = 100/i kms^^Mpc"^ flat 
ACDM cosmology, whilst all magnitudes are quoted in the Vega 
system. 



2 IMAGING 

2.1 Target fields 

The full VLT survey comprises 45 VIMOS pointings across nine 
quasar fields. In this paper we analyse an initial sample of 19 point- 
ings across 5 fields, where we have reduced and identified LBG 
spectra. The remaining LBG observations will be presented in a fu- 
ture paper. High-resolution optical spectra are available for all of 
the QSOs, which are at declinations appropriate for observations 
from the VLT at Cerro Paranal. The selected quasars for this paper 
are Q0042-2627 (z=3.29), SDSS J0124-^0044 (z=3.84), HE0940- 
1050 (z=3.05), SDSS J1201-H0116 (z=3.2 3) and PKS2126-158 
(z=3.28). Q0042-2627 has been observed bv lWilliger etaTI h996l) 
using the Argus multifibre spectrograph on the Blanco 4m tele- 
scope at Cerro Tololo Inter-American Observatory (CTIO) and as 
part of the Large B ri ght QSO Sur v ey (L BQS) using Keck/HIRES 
( iHewett et al1ll995h . Ipichon et al.l ( |2003|) observed HE0940-1050 
and PKS2126-158 using the Ultraviolet and Visual Echelle Spec- 
trograph ( UVES) on the VLT a nd SDSS J0124-I-0044 has been ob- 
served by IPeroux et all J2005l) also using UVES. Finally, SDSS 
J1201-I-0116 has been obse rved by the SDSS team using the 
SLOAN spectrograph and bv lO'Meara et alj ( l2007h using the Mag- 
ellan Inamori Kyocera Echelle (MIKE) high resolution spectro- 
graph on the Magellan 6.5m telescope at Las Campanas Observa- 
tory. 



2.2 Observations 

The imaging for our 5 selected fields was obtained using a com- 
bination of the MOSAIC Imager on the Mayall 4-m telescope at 
KPNO, the MOSAIC-II Imager on the Blanco 4-m at CTIO and 
VLT VIMOS in imaging mode. Q0042-2627, HE0940-1050 and 
PKS2126-158 were all observed at CTIO between January 2004 
and April 2005. JO 124-1-0044 and J 1201-1-01 16 were observed at 



KPNO in September 2001 and April 2006 respectively. All of these 
fields were observed with the broadband Johnson U (c6001) filter 
and the Harris B and R filters, except for JO 124-1-0044, which was 
observed with the Harris B, V and / broadband filters but not the 
Harris R. A full description of the observations is given in Table[T] 

We note that during the observations of the HE0940-1050 
field, there was a malfunction of one of the 8 CCDs leaving a gap 
of ~ 8' X 18' in the field of view. The remaining CCDs provided 
unaffected data however, which we use here. 

The MOSAIC Imagers each have a field of view of 36' x 36', 
covered by 8 2048 x 4092 CCDs. Adjacent chips are separated 
by a gap of up to 12" and we have therefore performed a dithered 
observing strategy for the acquisition of all our imaging data. For 
all observations we took bias frames, sky flats (durin g twilight pe- 
riods), dome flats and also observed iLandoltl <1992h standard-star 
fields with each filter on each night of observation for the calibra- 
tion process. 

In the Q0042-2627 and J1201-F01 16 fields, we also use imag- 
ing from the VLT VIMOS instrument with the broadband R filter. 
VIMOS consists of 4 CCDs each covering an area of 7' x 8', with 
gaps of 2' between adjacent chips. The fields were observed with 4 
separate pointings, with < 1' overlap between adjacent pointings. 

2.3 Data Reduction 

All data taken using the MOSAIC Imagers were reduced using the 
MSCRED package within IRAF, i n accordance with the NOAO Deep 
Wide-Field Survey guidelines of ljanuzzi et alj ( 120031) . Bias images 
were created using ZEROCOMBINE and dome and sky-flats were 
processed using CCDPROC. Removal of the "pupil-ghost" artifact 
was performed for the L'^-band calibration and science images using 
MSCPUPIL. 

The science images were processed using CCDPROC. Cosmic 
ray rejection was performed with CRAVERAGE in the early data- 
reductions (HE0940-1050 and PS2126-158), whilst in the later re- 
ductions, CRREJECT was used. The FIXPIX task was used to re- 
move marked bad-pixels and cosmic-rays from the images, using 
the interpolation setting. 

Deprojection of the images was performed using the MSCIM- 
AGE task, with optimization of the astrometry conducted using 
MSCCMATCH. Large-scale sky-variations were removed from sci- 
ence images using MSCSKYSUB and the resultant final images were 
combined using MSCIMATCH and MSCSTACK. 

For the HE0940-1050 and PKS 2126-158, short exposure 
imaging was obtained. These were used in the selection of QSO 
candidates (at brighter magnitudes than the LBG candidates) in 
these fields and were reduced and combined in the same way as 
the long exposure images described above. As there are typically 
only one or two short exposures per filter, the gaps between the 
CCDs still exist in the final short images, and no extra effort was 
made to remove blemishes by hand. 

The data reduction for the i?-band imaging from VLT VIMOS 
was performed using the VIMOS pipeline. Again bias frames were 
subtracted and the images were flat fielded using dome flats ac- 
quired on the night of observation. Individual exposur es were then 
depro jected and stacked using the SWARP software l IBertin et al.l 
l2002h . 



2.4 Pliotometry 

We performed object extraction using SEXTRACTOR, with a detec- 
tion threshold of 1 .2cr and a minimum object size of 5 pixels. Object 
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Table 1. Details of the imaging data acquired in each of our five target fields. Coordinates are given for the imaging centre, which is not necessarily the same 
as the position of the bright corresponding QSO. 



Field 


a 5 
(J2000) 


FaciUty 


Band 


Exp time 

(s) 


Seeing 


Depth 
50% comp. 


3fr 


Q0042-2627 


00:46:45 -25:42:35 


CT10/M0SA1C2 
VLTA'IMOS 


U 
B 
R 


12,600 
3,300 

235 


1.8" 
1.8" 
1.1" 


24.09 
25.15 
24.72 


26.16 
26.93 
25.79 



J0124+0044 


01:24:03 


+00:44:32 


KPNO/MOSAIC 


U 
B 
V 

I 


13,400 
2,800 
3,100 
7,500 


1.5" 
1.5" 
1.4" 
1.1" 



24.48 



25.60 
26.44 
26.14 

25.75 



HE0940-1050 


09:42:53 


-11:04:25 


CT10/M0SAIC2 


U 


29,000 


1.3" 


25.69 


26.75 










B 


4,800 


1.3" 


25.62 


26.66 










R 


2,250 


1.0" 


25.44 


26.24 


J1201+0116 


12:01:43 


+01:16:05 


KPNO/MOSAIC 


U 


9,900 


1.6" 


24.50 


26.11 










B 


6,000 


2.4" 


24.43 


26.56 








VLTA'IMOS 


R 


235 


0.7" 


25.47 


26.24 


PKS2126-158 


21:29:12 


-15:38:42 


CT10/M0SA1C2 


U 


26,400 


1.3" 


25.08 


26.97 










B 


7,800 


1.6" 


24.94 


27.49 










R 


6,400 


1.5" 


24.65 


26.79 



detection was performed on the 7?-band images and fluxes were 
calculated in all bands using Kron, fixed-width (with a diameter 
of twice the image seeing FWHM) and isophotal width apertures. 
Zeropoints for each of the observations were calculated from the 
Landolt standard-star field observations made during the observing 
runs and we correct the photometry for galactic extinction using the 
dust maps of ISchlegel et alj ( 119981) . Each of the standard-star field 
images were processed using the same method as for the science 
frames. The depths reached in the U, B and R bands for each field 
are given in table [T] We quote the 3cr depths, which give the limit 
for detecting an object 5 pixels in size with a signal of 3 x the back- 
ground RMS detection, and the 50% completeness level. The 50% 
completeness levels are calculated by systematically placing sim- 
ulated point-source objects in the final stacked images at different 
magnitudes. The 50% level is then the magnitude at which we are 
able to recover 50% of simulated sources. 

The U, B and R number counts from the 4 fields are plotted 
in Figs.[T]to[3] In general the counts turnover at ~ 0.5mag brighter 
than the 50% completeness limits, consistent with the counts being 
dominated by extended sources (whilst the completeness limits are 
estimated using simul ated point-sources). W e plot for comparison 
the number counts of iMetcalfe et alj ( 1200 ID . All counts are from 
our MOSAIC data except for the R band counts of Q0042-2627 
and JI201+01I6, which are from the VLT VIMOS. The imaging 
in the JI201+01I6 field was taken during relatively poor seeing 
conditions during observations at CTIO and so reaches shallower 
depths than the other fields. For these plots, stars have been re- 
moved using the SEXTRACTOR CLASS.STAR estimator with a limit 
of CLASS.STAR < 0.8. 



2.5 Selection Criteria 



10= 



We perform a photometric selection based on that of ISteidel et al.l 
( Il996ll2003l) . but applied to the U, B and R band ima ging avail- 
able from our imaging survey. As in lsteidel et al.l J2003h our selec- 
tion takes advantage of the Lyman-Break at 912A and the Lya- 
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Figure 1. C-band number counts from the four fields Q0042-2627 (black 
crosses), HE0940-1050 (dia monds), J1201+0116 (t riangles) and PKS2126- 
158 (squares). The counts of lMetcalfe eTal] J200lh from the William Her- 
schel Deep Field are shown for comparison (red crosses). 



forest passing through the [/-band and into the B-band in the 
redshift range 2.0 < z < 3.5. To est ablish the selection in the 
Vega UBR system, we convert from the lSteidel et alJ(l2003h selec- 
tions u sing the photometric transformations of lSteidel & Hamilton! 
( 119930 . moving from the UnGTZ AB system to the Johnson- 
Morgan/Kr on-Cousins Vega photome try. The approximate trans- 
formations dSteidel & Hamiltoril 19931) are as follows: [/„ = t/ + 
0.75, G = B - 0.17 and 7^ = i? + 0.14 and transform the 
ISteidel et alj ( l2003h selection to {B - R) s^ 1.51 and {U - B) ^f 
(B-R)- 0.23. 
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Figure 2. _B-band number counts from the four fields Q0042-2627 (black 
crosses), HE0940-1050 (dia monds). J1201+0116 (t riangles) and PKS2126- 
158 (squares). The counts of lMetcalfe et al.l J200lh from the Wilham Her- 
schel Deep Field aie shown for comparison (red crosses). 
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Figure 3. ij-band number counts from the four fields Q0042-2627 (black 
crosses). HE0940-1050 (dia monds). J1201+0116 (t riangles) and PKS2126- 
158 (squares). The counts of lMetcalfe et al.l J200ll) from the Wilham Her- 
schel Deep Field ai'e shown for comparison (red crosses). 



We also tak e into account model colo ur tracks calculated us- 
ing GALAXEV JBmzual & CharlollbOOJh . The tracks are shown 
in Figs. |4] and |5] (solid black curves). We use a Salpeter initial 
mass-function, assuming solar metallicity with a galaxy formed at 
z = 6.2 (i.e. with an age of 12.6 Gyr at z = 0) and a r = 9 Gyr 
exponential SFR. The three different curves show the effect of dust 
extinction with a model given by (left to right) r^ = 0.5, r^ = 1.0 
and Tv = 2.0, where r^ = 2.0, where r^ is the effective absorption 



JCharlot & Fall[ 20000. The models agree well with the transfor- 
mation of the ISteidel et alj ( I2003I) selection criteria, although the 
dustier models do suggest a greater extensio n of the z > 3 pop- 
ulation to higher values of (B - R) than the lSteidel et alj ( |2003|) 

criteria. 

Based on the models and the lSteidel et alj ( 120031) criteria, we 
develop a number of selection criteria in the UBR system. The key 
modification s that we make from our initial colour-cut estimates 
based on the ISteidel et alj ( 120030 cuts are to extend the selection 
further redwards in [B — R) and to align the {U — B) — {B — R) 
axis with the stellar locus in the UBR plane, which has a slope of 
((7 - B) ~ 1.25{B - R). We note that the first of these modi- 
fications ri sks increasing the number of contaminants in the form 
of M-stars ISteidel & Hamiltonll 19930 and the second increases the 
risk of contaminants in the form of lower redshift galaxies. How- 
ever, given the large number of slits available to us with the VLT 
VIMOS spectrograph, we deem the risk of increased levels of con- 
tamination acceptable, whilst extending the colour-cuts can allow 
the observation of dusty z > 3 objects as well as 2 « 3 galax- 
ies which may be scattered out of the primary selection area due 
to photometric errors on these faint objects. As such we use four 
selection criteria with different priorities for spectroscopic obser- 
vation (taking advantage of the object priority system in arranging 
the VIMOS slit masks). These selection criteria are as follows: 

• LBGJ'RIl 

(i) 23 < fi < 25.5 

(ii) U - B >0.5 

(iii) B - R< 0.8(1/ - B)+ 0.6 

(iv) B- R<2.2 

• LEGJ^RE 



(i) 23 < R< 25.5 
(ii) U - B >0.0 
(iii) B - R< 0.8{U - 
(iv) B- R< 2.8 

• LBG^RB 



B) + 0.8 



(i) 23 < H < 25.5 

(ii) -0.5 <U ~ B <0.0 

(iii) B ~ R< 0.8(1/ ~ B) + 0.6 

• LBGJDROP 

(i) 23 < H < 25.5 
(ii) No U detection 
(iii) B- R<2.2 

LBG_PRI1 is our primary sample and selects candidates that 
are expected to be the most likely 2.5 < 2 < 3.0 galaxies. The 
LBG_PRI2 sample targets objects with colours closer to the main 
sequence of low-redshift galaxies than the LBG_PRI1 objects. This 
sample is therefore expected to include a greater level of con- 
tamination from low redshift galaxies. In addition, based on the 
path of the evolution tracks in Figs. |4] and \5\ we also expect the 
z > 2.5 population that this selection samples to have, on aver- 
age, a lower redshift than the LBG_PRI1 sample. The next selec- 
tion sample, LBG_PRI3, takes this further and is intended to target 
a 2.0 < z < 3.0 galaxy redshift based on the evolution tracks. Fi- 
nally, we select a sample of [/-dropout objects (LBGJDROP) with 
detections in only our B and R band data. 

In none of the above samples do we attempt to remove stellar- 
like objects due to the risk of losing good LBG candidates. The 
half-light radius of 2 ~ 3 LBGs has been shown to be on average 
0.4" and so will not be resolved in our data, which is mostly taken 
under conditions of > 0.8" seeing. 
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We apply these selection criteria to four of our QSO fields: 
Q0042-2627, HE0940-1050, J1201+0116 and PKS2126-158. The 
candidate selection for the JO 124+0044 field was performed sepa- 
rately and is discussed in Bouche & Lowenthal (2004). Figs.|4]and 
|5]show the four selection criteria applied to these four fields. The 
selection boundaries are shown by the red, green and blue lines for 
the LBG_PRI1, LBG_PRI2 and LBG_PRI3 selections respectively. 
Objects selected as candidates by each criteria set are shown by 
red, green, blue and cyan points for the LBG_PRI1, LBG_PRI2, 
LBG_PRI3 and LBG_DROP selections respectively. The grey con- 
tours in each plot show the extent of the complete galaxy population 
in each of the fields. 

Returning to the depths of our fields, we now compare these 
to t hose of previous st udies in the selection of LBGs. We note 
that lSteidel et al.l ( l2003h used photometry with mean Icr depths of 
((t([/„)) = 28.3, (a(G)) = 28.6 and (o-(7^)) = 28.0, whilst their 
imposed IZ band limit was TZ = 25.5. Using t he trgin sformations of 
ISteidel & HamiltonI h993h , the lSteidel et"S ( l2003h la limits cor- 
respond to [/ = 27.55, B = 28.77 and R = 27.86 in the Vega 
system. Comparing this to the average depths in our own fields, we 
have mean 3cr depths of [/ = 26.2, B = 26.8 and R = 26.3, 
which equate to la depth s of t/ = 27.4, B = 28.0 and R = 27.5, 
largely comparable to the lSteidel et al.l <2003h imaging data. 

The numbers of objects selected by each selection for each 
field are given in Table [2] These candidate selections were used 
as the basis for the spectroscopic work which is described in the 
following sections. 



2.6 QSO Candidate Selection 

At redshifts of 2: ~ 3, the observed optical spectra of QSOs and 
galaxies exhibit similar shapes, both being heavily influenced by 
the Lyman break feature. We therefore add to our targets a number 
of QSO candidates in each field (except J0124+0044) using the fol- 
lowing selection, which is closely based on our high-priority 2 ~ 3 
LBG selection: 

(i) CLASS.STAR > 0.8 

(ii) (7 - B > 0.5 

(iii) B - R< 0.8((7 ~ B) + 0.8 

(iv) 0.0 < B - i?, < 2.2 

The magnitude limits used with this selection were 20 < i? < 
23 in the Q0042-2627 and J1201+0116 fields and 18 < i? < 22 in 
the HE0940-1050 and PKS2126-158 fields for which we had ob- 
tained shallow imaging and could therefore select brighter objects 
more reliably. 

As with the LBGs, QSOs at z > 2 may be sele cted by the pas- 
sage of the Lyman-break through the [/-band (e.g. [Richards et al.l 
120091) . This selection is therefore based on the LBG selection, but 
constrained to brighter magnitudes and stellar-like objects only. 
This selection gives 71, 39, 15 and 38 QSO candidates in the 
Q0042-2627, HE0940-1050, J1201+0116 and PKS2126-158 fields 
respectively. Note that only a small number of these have actually 
been observed spectroscopically as the LBG candidates remained 
the higher priority. 



3 SPECTROSCOPY 

3.1 Observations 

We observed our LBG candidates using the VIMOS instrument on 
the VLT UT3 (Melipal) between September 2005 and March 2007. 



As described earlier, the VIMOS camera consists of four CCDs, 
each with a field of view of 7' x 8', arranged in a square config- 
uration, with 2' gaps between the field-of-views of adjacent chips. 
Each observation therefore covers a field of view of 16' x 18' with 
224 arcmin^ being covered by the CCDs. The instrument was set 
up with the low-resolution blue grating (LR_Blue) in conjunction 
with the OS_Blue filter, giving a wavelength coverage of 3700A to 
6700A and a resolution of 180 with 1" slits, corresponding to 28A 
FWHM at 5OOOA. The dispersion with this setting is 5.3A per pixel. 
We note that this configuration also projects the zeroth diffraction 
order onto the CCDs. 

Given the size of our imaging fields (36' x 36') it was possible 
to target 4 distinct sub-fields with the VIMOS field of view. We 
have therefore observed a total of 19 sub-fields across our 5 fields, 
i.e. 4 sub-fields in each field except for HE0940-1050 in which only 
3 sub-fields were achievable due to the CCD malfunction during 
the imaging observations. Each sub-field was observed with 10 x 
1, 000s exposures, apart from sub-field three of the PKS2126-158, 
which was observed with only 4x1, 000s due to time constraints 
in the VIMOS schedule. All observations were performed during 
dark time, with < 0.8" seeing and < 1.3 air mass. 

Slit masks for each quadrant of each sub-field were designed 
using the standard VIMOS mask software, VMMPS. We used min- 
imum slit lengths of 8", which equates to 40 pixels given the pixel 
scale of 0.205"/pixel. With the effectively point-like nature of our 
sources and our maximum seeing constraint of 0.8" this allows us 
a minimum of ~ 7" for sky spectra per slit (with which to per- 
form the sky-subtraction when extracting the spectra). Using the 
VMMPS software with the LR_Blue grism we were able to tar- 
get up to ~ 60 — 70 objects per quadrant (i.e. « 250 objects 
per sub-field), depending on the sky density of the candidate ob- 
jects. For the spectroscopic observations, we predominantly used 
the selections as given in section [23] however to optimize the spec- 
troscopic observations some flexibility was employed in including 
small numbers of objects outside the selection criteria. However, 
we note that the LBG_PRI3 selection was not employed in the spec- 
troscopic observations in the first observations (i.e. the observations 
of HE0940-1050 and PKS2126-158), whilst the magnitude limit 
used for selecting objects to observe for later fields was reduced 
from R = 25.5 Xo R = 25. The total number of spectroscopically 
observed objects was 3,562. 



3.2 Data reduction 

Bias frames were obtained by the VLT service observers at the be- 
ginning of each night of observations. Lamp-flats were also taken 
with each of the masks with the observation setup in place (i.e. the 
OS_Blue filter and LR_Blue grism). These were also taken by the 
service observers at the beginning of each night's observation. Arc 
frames were taken during the night with each of the masks with the 
LR_Blue grism and OS_Blue filter. 

Data reduction was performed using the VIMOS pipeline soft- 
ware, ESOREX. Firstly the bias frames were combined to form a 
master bias using VMBIAS. The flat frames were then processed 
and combined using the VMSPFLAT recipe. VMSPCALDISP was 
then used to process (bias subtract and flat-field) the arc lamp expo- 
sure and to determine the spectral distortions of the instrument. We 
measured a mean RMS on the inverse dispersion solution (IDS) of 

2.3 ± 0.6 A. With the bias, flat and arc exposures all processed, the 
object frames were reduced and combined using the VMMOSOB- 
SSTARE recipe to produce the reduced 2-D spectra. The spectra 
have not been fully flux calibrated, however we have applied the 
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Figure 4. Our selection criteria in UBR colour space shown for the Q0042-2627 (left) and HE0940-1050 (right). The red line and points show the LBGJ'RIl 
selection, the green line and points show the LBGJ'RIl selection, the blue line and points show the LBG_PRI3 selection and the cyan line at (7 — _B = 4.5 
shows the LBG_DROP selection. The grey contours show the entire galaxy population in the fields. The black lines show the galaxy evolution model for a 
galaxy with a r = 9Gyr exponential SFR formed at z = 6.2 and are labelled with values of observed redshift from z = 3.83 to z = 0. 

Table 2. Number of candid ate high redshift objects in each of the selected fields. Note that candidates in the J0124+0044 were selected as described in 
iBouche & Lowenthal <2004f) and not using the four selection criteria sets described in this paper 



Field 


LBGJ'RIl 


LBG.PRI2 


LBGJ'RI3 


LBGJDROP 


Total 


Q0042-2627 


1,366 


1,381 


650 


1,390 


4,787 


J0124+0044 










3,679 


HE0940-1050 


1,646 


2,249 


741 


1,042 


5,678 


J1201+0116 


477 


487 


469 


606 


2,029 


PKS2126-158 


1,380 


2,119 


713 


667 


4,879 


Total 


4,869 


6,236 


2,573 


3,705 


21,062 


Observed spectroscopically 


730 


569 


256 


999 


2,554 



master response curves for the LRJilue grism to correct for the 
effects of the grism as a function of wavelength. 

We extract the 1-D spectra using purpose- written IDL rou- 
tines. For each spectrum, we first fit the shape of the spectrum 
across the slit. This is implemented by binning the 2-D aperture 
along the dispersion axis and then fitting a Gaussian profile to each 
bin to find the centre of the object signal in each bin. We then fit 
the resultant spread in the central pixel with a 4th order polynomial 
function. We then lay an object aperture with a width of Hap pixels 
over the object and a sky aperture covering all of the usable sky re- 
gion in the slit. The object and sky spectra are then taken as being 
the mean over the widths of their respective apertures. Finally, we 
subtract the sky spectrum from the object spectrum to produce the 
final object spectrum. The dominant remaining sky-contamination 



after sky-subtraction were the strong sky emission lines [OI]5577 
A [Nal]5890 A and [OIJ6300 A. 

We estimate the signal-to-noise by taking the RMS of the sky 
aperture in each wavelength bin and dividing by y/Uap, where Uap 
is the width of the aperture used to extract the 1-D spectrum of a 
given object. Fig. [6] shows the mean signal-to-noise per resolution 
element (i.e. 28A) in the wavelength range 4100A< A <5300A 
in our sky- subtracted spectra as a function of source 7?-band mag- 
nitude. The selected range covers many of the key emission and 
absorption lines exhibited in LBGs in the redshift range 2.5 < z < 
3.5, whilst excluding the strong sky lines. The points in Fig.|6]show 
the mean spectrum SNR per resolution element, whilst the error 
bars show the standard deviation within each bin. In the faintest 
bin (25.25 < R < 25.5), we achieve a mean continuum signal-to- 
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Figure 5. As in Fig.|4]but for the J1201+01 16 and PKS2126-158 fields (left to right). 



Table 3. Details of the spectroscopic data acquired in each of our five target fields. Coordinates are given for the targeting centre of each sub-field. 



Field 


Sub-field 


o 


<5 


Dates 


Exp time 


Seeing 






(J2000) 


(J2000) 




(s) 






Q0042-2627 


fl 


00:45:11.14 


-26:04:22.0 


8-10,15/08/2007 


10, 000 


0.6- 


- 1.0" 


Q0042-2627 


f2 


00:43:57.30 


-26:04:22.0 


18-19/08/2007 & 5-6/09/2007 


10, 000 


0.9- 


^ 1.0" 


Q0042-2627 


f3 


00:45:10.35 


-26:19:06.9 


11-12/09/2007 


10, 000 


0.9- 


- 1.0" 


Q0042-2627 


f4 


00:43:55.97 


-26:19:16.1 


7,10/09/2007 


10, 000 


0.9- 


- 1.0" 


JO 124+0044 


fl 


01:24:41.82 


+00:52:18.8 


1-2,4/11/2005 


10, 000 


0.8- 


- 0.9" 


JO 124+0044 


f2 


01:23:32.06 


+00:52:13.1 


5,29,31/10/2005 


10, 000 


0.6- 


- 1.0" 


JO 124+0044 


f3 


01:23:31.29 


+00:37:02.0 


19-20/09/2007 


10,000 


0.8- 


- 1.0" 


JO 124+0044 


f4 


01:24:41.86 


+00:36:51.4 


4/12/2005 & 22/08/2006 


10,000 


0.8- 


- 0.9" 


HE0940-1050 


fl 


09:42:08.02 


-11:08:14.2 


26-27,29/01/2006 


10,000 


0.5- 


- 0.8" 


HE0940-1050 


f2 


09:43:21.53 


-11:08:35.0 


30-31/01/2006, 1,25/02/2006 & 1/03/2006 


10, 000 


0.5- 


^ 1.0" 


HE0940-1050 


f3 


09:43:21.58 


-10:54:31.8 


14,19/12/2007 & 31/01/2008 


10, 000 


0.6- 


- 1.0" 


JI201+0116 


fl 


12:02:14.01 


+01:09:09.9 


13-15/04/2007 & 17/04/2007 


10, 000 


0.6- 


- 1.0" 


J1201+0116 


12 


12:01:10.01 


+01:09:09.9 


23/04/2007 & 8,11,14/05/2007 


10, 000 


0.4- 


- 0.9" 


J1201+0116 


f3 


12:01:10.04 


+01:24:09.8 


16-17/05/2007 


10, 000 


0.5- 


- 0.9" 


J1201+0116 


f4 


12:02:14.07 


+01:24:08.0 


18/05/2007 & 6,8,10/02/2008 


10, 000 


0.6- 


-0.7" 


PKS2126-158 


fl 


21:29:59.57 


-15:31:30.2 


17/08/2006 & 1,21-26/09/2006 


10, 000 


0.7- 


- 1.0" 


PKS2126-158 


f2 


21:28:46.20 


-15:31:29.9 


5-6/08/2005 


10, 000 


0.6- 


- 1.0" 


PKS2126-158 


f3 


21:30:00.41 


-15:47:18.3 


27/09/2006 


4,000 


0.8- 


- 1.0" 


PKS2126-158 


f4 


21:28:46.27 


-15:47:11.9 


9-11,25,29/08/2005 


10,000 


0.7- 


- 0.9" 



noise of ~ 3.5. This rises to a continuum signal-to-noise 
our brightest objects (23 < i? < 23.25). 



9 for 3.3 Object Identification 

We perform the object identification for each sht individually by 
eye. Given the wavelength range covered by the LR_Blue grism 
combined with the redshift range of our targets, 2 < z < 3.5, there 
are several key spectral features that facihtate the identification of 
those targets. These are primarily: 
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Figure 6. Mean signal-to-noise per resolution element (28A) in the wave- 
length range 4100A< A <5300A as a function of i?-band magnitude in 
our VLT VIMOS spectra with integration times of 10,000s. 



• Lyman limit, 91 2A; 

• Ly/3 emission/absorption, 1026A 

• OVI 1032A, 1038A; 

• Lya forest, <1215.67A; 

• Lya emission/absorption, 1215. 67A; 

• Inter-stellar medium (ISM) absorption lines: 

- Sill 1260.4A; 

- Ol+Sill 1303 A; 

- CII 1334A; 

- SilV doublet 1393 A & 1403 A; 

- Sill 1527A; 

- Fell 1608A; 

- AlII 1670A; 

• CIV doublet absorption/emission, 1548-1550A. 

The most prominent of these features is most frequently the 
Ly a emission/ab s orptio n feature at 1215A. However, as discussed 
bv lShaplev et alj ( 120031) , the observed optical (rest-frame UV) ab- 
sorption and emission features are thought to originate from an out- 
flowing shell of material surrounding the core nebular region of the 
galaxy. These features do not therefore represent the redshift of the 
rest-frame of the galaxy but in fact of these outflows. 

For each confirmed LBG we measure independently the red- 
shift of the Lya emission/absorption feature and the redshift of the 
ISM absorption features. In order to measure the Lya redshift, we 
fit the feature with a Gaussian function allowing the amplitude, 
central wavelength and width to be free parameters. From these 
we determine the redshift and line-width of the feature. We note 
that absorption blue- wards of the emission wavelength produces an 
asymmetry in the observed emission line, however given the mod- 
est resolution of our observations the Gaussian fit is preferred to 
any more complex asymmetric fitting to the emission line. 

We have performed an estimate of the accuracy of our red- 
shift results by repeating the spectral line fitting method with mock 
spectra. Each mock spectrum consists of a single Gaussian emis- 
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Figure 7. Estimate of the accuracy of the Gaussian line-fitting based on 
iteratively fitting mock spectra with Gaussian random noise. The open cir- 
cles show the results of applying the fitting method to a single emission 
line spectrum with a range of signal-to-noise (where the signal-to-noise is 
defined as the ratio between the peak signal and the width of the Gaussian 
noise). The blue triangles show the result of the same method applied to 
a simple absorption line spectrum including the ISM lines: Sill (1260A), 
Ol-l-Sill (1303A), CII (1336A) and SilV (1393A, 1402A). 



sion line (i.e. / = Ae~'^^~^°'> 1"^" ) at a random redshift in the 
range 2.5 < z < 3.5 and a FWHM of 1680 kms"^ corresponding 
to a Gaussian width of a = 850 kms"^ (equivalent to the resolu- 
tion of the instrument). Gaussian random noise was then added to 
the basic emission line shape to give the required signal-to-noise. 
For each mock spectrum, we then performed the Gaussian fitting, 
iteratively performing the process for a total of 10* mock spectra at 
a given signal-to-noise. The difference between the input redshift 
and the Gaussian line fitting redshift was then measured for each of 
the iterations and the error estimated from the distribution of this 
difference in input and measurement. The process was repeated, 
increasing the emission line peak flux from 1 to 20x the Gaussian 
noise width. 

The results are given in Fig.|7] where the measured accuracy is 
plotted as a function of the calculated signal-to-noise (red circles). 
Further to this, we measure the distribution of Lya emission peak 
signal-to-noise in our galaxy sample, which is shown in Fig.[8]as a 
percentage of the total number of LBGs exhibiting Lya emission. 
If we now compare these two plots, we see that ~ 90% of our emis- 
sion line LBGs have an emission line signal-to-noise of > 3, which 
suggests that 90% of the Lya emission line redshifts have velocity 
errors of less than ~ 550 kms^^. Further, the median Lya emis- 
sion line signal-to-noise is ~ 5.5 which gives a velocity error of 
~ 400 kms~^. Our higher quality spectra (i.e. the top 20%) how- 
ever, are estimated to achieve velocity errors on the Lya emission 
line redshifts as small as ~ 200 kms~^. 

Where feasible, we also attempt to measure the redshift of the 
ISM absorption lines based on the Sill, Ol-l-Sill, CII and SilV dou- 
blet (despite being a mixture of high and low ionization lines we 
note t hat they are all mea sured to have comparable velocity off- 
sets in lShapIev et alj2003l . at least within the resolution constraints 
afforded by our observations). We primarily use absorption lines 
between 1215A < A^est ^ 1500A as these remain within the 
wavelength coverage of the low-resolution blue grism over the full 
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Figure 8. The distribution of Lya emission line (red stars) and ISM absorp- 
tion line (blue circles) signal-to-noise measurements in our LBG sample. 
The calculated signal-to-noise is the ratio between the emission/absorption 
line peak (after subtracting the continuum) and the measured noise. The fi- 
nal ISM signal-to-noise value is taken as the median of the calculated values 
for the ISM Unes used. See Fig.|2]for the estimated velocity errors based on 
the feature signal-to-noise. 



redshift range (i.e. 2 < z < 3.5) of our survey. Measuring the indi- 
vidual absorption lines in most of our spectra is difficult given the 
SNR of the absorption features in our spectra, however our ability 
to estimate the redshift of the ISM lines can be greatly improved by 
attempting to determine the mean ISM redshift by fitting the five 
lines simultaneously. 

To evaluate this method we repeat the iterative error analysis 
performed for the Lya emission line fitting, but fitting five absorp- 
tion lines (with aisM = 850 kms^^) simultaneously. Again we 
measure the offset between the input redshift and the output redshift 
measured from the Gaussian line fitting. The result is again plotted 
in Fig. [7] (blue triangles), whilst the distribution of ISM signal-to- 
noise measurements in the data is again given in Fig. [8] This sug- 
gests that we may reasonably expect a significant improvement in 
the estimated redshift compared to measuring just a single line. We 
now predict an accuracy of « 200 kms~^ at a signal-to-noise of 
~ 3, which based on Fig.[8]accounts for 55% of our sample. 

With the Lya and ISM redshifts determined, we estimated the 
intrinsic redshifts, Zjnt, of our LBG sample using the relations of 
lAdelberger et aLlfeOOSah . These relations were derived from a sam- 
ple of 138 LBGs observed spectroscopically in both the optical and 
the near infrared and are based on the offsets found between the 
Lya plus ISM lines and the nebular emission lines, [OII]3727A, 
H/?, [OIII]5007A and Ha. These lines are all associated with the 
central star-forming regions of LBGs as opposed to the outflowing 
material and are thus expected to be more repres entative of the in- 
trinsic redshift of a given LBG. The relations of lAdelberger et al.l 
( l2005ar) that we use here are as follows: 

For LBGs with only a redshift from the Lya emission line we 
used: 



2int = 2LyQ - 0.0033 - 0.0050(2Lya " 2.7) 



(1) 



For objects with Lya absorption and a measurement of zism we 
used: 



2int 



2ISM + 0.0022 + 0.0015(2isM - 2.7) 



(2) 



And for objects with redshifts measured from both the Lya emis- 
sion line and the ISM absorption lines we used: 

2int = 5 + 0.070A2 - 0.0017 - 0.0010(5 - 2.7) (3) 

where z is the mean of the Lya redshift (zLya) an d the ISM absorp- 
tion lin e redshift (zism) and Az = ZLya — ztsM. lAdelberger et al.l 
( l2005d) quote rms scatters of g^ = 0.0027 (200 kms"^), 0.0033 
(250 kms~^) and 0.0024 (180 kms"^) respectively for each of the 
above relations based on their application to their optical and IR 
spectroscopic sample of LBGs. 

As well as 2 « 3 galaxies, our selection also samples a num- 
ber of contaminating objects. These consist of low-redshift emis- 
sion line galaxies (identified by [OII]3727A, H/3, [OIII]5007A and 
Ha emission), low-redshift Luminous Red Galaxies (LRGs - iden- 
tified by [OIIJ3727A emission, Ca H, K absorption and the 4000A 
break) and faint red stars (mostly M and K-type stars). We show 
examples of the spectra of several LBGs and contaminant low- 
redshift galaxies taken with the VLT VIMOS in this survey in Fig.|9| 
(note that these are not flux-calibrated spectra). 

All identified objects, including stars and low-redshift galax- 
ies, were assigned a quality rating, q, based on the confidence of 
the identification. The value of q was assigned on a scale of to 1, 
with 1 being the most confident and being unidentified. All ob- 
jects with q < 0.5 were rejected as spurious identifications and are 
not included in the spectroscopic catalogue used in the remainder 
of this work. LBGs were generally classified as follows: 

• 0.5 - Lya emission or absorption line evident plus some 
'noisy' ISM absorption features. 

• 0.6 - Lya emission or absorption plus some ISM absorption 
features. 

• 0.7 - Lya emission or absorption plus most ISM absorption 
features. 

• 0.8 - Clear Lya emission or absorption plus all ISM absorption 
features. 

• 0.9 - Clear Lya emission or absorption plus high signal-to- 
noise ISM features. 

With this classification scheme, we have identified 392, 254, 
170, 111 and 93 2 > 2 galaxies with q =0.5, 0.6, 0.7, 0.8 and 0.9 
respectively. 



3.4 Sky Density, Completeness & Distribution 

We summarize the numbers of objects observed in Table |4| 
Our mean sky density for successfully identified LBGs is 
0.24arcmin~^, whilst the percentage of z > 2 galaxies in the en- 
tire observed sample (the success rate given in table |4]l is 27.5%. 
The remaining observed objects are a mix of low-redshift galaxies, 
stars and unidentified objects (generally very low-signal to noise 
spectra). In the worst case field (J1201-I-0116), we have a greater 
number of low-redshift galaxies than high redshift detections. We 
attribute this to the relatively poor depth of the imaging observa- 
tions in this field. We also note that the PKS2126-158 field is at 
a relatively low galactic latitude and thus was a higher proportion 
of contamination by galactic stars. However, the field still shows a 
high proportion of 2; > 2 galaxies. 

In Fig.[TO|and Table|5]we summarize the redshift distributions 
of each of our sample selections in our observed fields. The overall 
redshift distribution across all fields is shown in the bottom panel of 
Fig. [To] with the black histogram showing the redshift distribution 
from IJBVI selected objects from J0124-I-0044 and the red, green, 
blue and cyan histograms showing the LBG_PRI1, LBG_PRI2, 
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Figure 9. Example spectra taken using 10,000s integration time witli the LR_Blue grism on tiie VLT VIMOS instrument. The top two spectra are examples 
of contaminating low-redshift galaxies. The remaining 12 panels show LBG spectra exhibiting both Lya emission and absorption over the redshift range 
2 < 2 < 3.5. ISM lines are also clearly identifiable in the individual LBG spectra as is the Lyman limit. Both galaxy redshift and apparent i?-band magnitude 
(Vega) are quoted for each object. Note that all the above spectra have been binned to Si 16A. 

Table 4. Summary of objects identified in the VLT VIMOS observations. The success rate is the number of successfully identified LBGs divided the total 
number objects observed. Example spectra of the high-redshift and low-redshift galaxies are shown in Fig.|9] All 10 identified z > 2 QSO spectra are provided 
in Fig. 1161 



Field 


Subfields 


Slits 


Galaxies 
z > 2 


QSOs 

z> 2.0 


Galaxies 
z< 2.0 


Stars 


Success rate 


Q0042-2627 


4 


876 


264 (0.29arcmin-2) 


1 


106 


5 


30.1% 


10124+0044 


4 


832 


264 (0.29arcmin~2) 





54 


18 


31.7% 


HE0940-1050 


3 


501 


169 (0.25arcmin~^) 


1 


48 


36 


33.7% 


J1201+0116 


4 


699 


120(0.13arcmin-2) 


5 


144 


72 


17.2% 


PKS2126-158 


4 


654 


203 (0.23arcmin-2) 


3 


49 


126 


31.0% 


Total 


19 


3562 


1020(0.24arcmin-2) 


10 


401 


257 


28.6% 
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Figure 10. Differential redsliift distribution in each of our fields and 
summed over all fields. We show the number counts split by selec- 
tion criteria: LBGJDROP (cyan histograms), LBG_PRI1 (red histograms), 
LBG_PRI2 (green histograms) and LBG_PRI3 (blue histograms). The mean 
redshifts for each selection are given in table [5] 



LBG_PRI3 and LBG_DROP selections respectively. The overall 
mean redshift for our confirmed LBG sample is 5 = 2.85 ± 0.34. 
It is evident from the redshift distributions that the separate selec- 
tion sets give slightly differing (but overlapping) segments in red- 
shift space. As may be expected, the LBGJDROP selection is the 
most biased towards the higher end of our redshift distribution, with 
an overall mean redshift across all our samples of 5 = 2.99. The 
LBG_PRI1 selection provides a redshift range of 2.90±0.32, whilst 
the LBG_PRI2 and LBGJRB give comparable redshift distribu- 
tions of 2.67 ± 0.26 and 2.67 ± 0.31 respectively. We also show 
the redshift distributions for each individual field in the top five 
panels of Fig. [H with the LBG_PRI1, LBG_PRI2, LBGJRD and 
LBGJDROP identically to that in the 'all fields' plot. In each field 
we again see that the LBGJRIS and LBGJ'Rn selections pref- 
erentially select the lowest redshift ranges followed by LBGJRII 
and LBGJDROP showing the highest redshift range (although this 
is less pronounced in the J 1201-1-01 16 field in which the imaging 
depths were least faint). 

We illustrate the distribution of our spectroscopic LBG sample 
in each of our 5 fields in Fig.[TT] The fields are ordered by R.A. top 



to bottom and all identified z > 2 galaxies (filled blue circles) 
are shown along with all known z > 2 QSOs identified from the 
NASA Extragalactic Database. We also plot the positions of QSOs 
identified in our VIMOS obs ervations and AAOme ga QSO survey, 
which is described further in lCrighton et al.l OOlOl) . 

In Fig. [12] we plot the number of identified LBGs in magni- 
tude bins for each of our fields. The filled histograms show the 
cumulative numbers of successfully identified objects (including 
interlopers as well as 2 > 2 galaxies) split by their selection 
criteria. LBGJDROP selected objects are shown by the cyan his- 
togram, LBGJ'RIl by the red histogram, LBGJ'Rn by the green 
histogram and LBGJ'RD by the blue histogram. The distribution 
of all spectroscopically observed objects is given by the solid line 
histogram in each case. As the J0I24-I-0044 objects were not se- 
lected using the same selection criteria, these are simply left as a 
single group shown by the filled black histogram. In all fields, we 
see that we are successfully identifying objects down to the magni- 
tude limit of R= 25.5 (1= 25 in the case of JO 1 24-1-0044), although 
a significant number of objects remain unidentified in each field at 
the fainter magnitudes as spectral features become more difficult to 
discern in the spectra. We note also that the shapes of the overall 
magnitude distributions are biased more towards brighter objects 
in the Q0042-2627 and J 1201-1-01 16 fields in which a greater num- 
ber of LBGJ'RD objects are included (and also the imaging depths 
achieved in these fields are shallower than in the other fields). 

In Fig. [T3] we show the number counts of our photometri- 
cally selected LBGs (open red circles) and the estimated number 
counts of LBGs (filled red circles) derived from the candidate num- 
ber counts and the success rate as a function of magnitude (i.e. 
the number of confirmed LBGs divided by the number of observed 
candidates). At faint magnitudes we correct the counts for incom- 
pleteness in the spectroscopic observations, however we have not 
made any correction f or incompleteness in the original photometry. 
The number counts of ISteidel et al.l 1 120030 are also plotted, show- 
ing their candidate number counts (open blue triangles) and number 
counts corrected for contamination (filled blue triangles). The two 
data-sets show good agreement over the magnitude ranges sampled. 



3.5 Velocity Offsets and Composite spectra 

The galaxy spectra contain a wea lth of information as illustrated by 
the work of lShaplev et al.l ( 120031) . We now look at how our spectra 
compare to previous work in terms of the velocity offsets between 
the different spectral features. For the galaxies that exhibit both 
measurable Lya emission and ISM absorption lines, we calculate 
the velocity offsets between these lines, A« = Vem ~ Vabs- The 
distribution of Aii for our galaxy sample is shown in Fig. [14] The 
distribution of velocity offsets exhibits a strong peak with a mean 
of (An) = 625 with a disp ersion of 510 k nis~^ . This compares to 
a value measured bv .Shaplev et al.l J2003h of 650 kms"^ . 

We have produced composite spectra in several Lya equiva- 
lent width bins in order to produce spectra with increased signal- 
to-noise compared to the individual galaxy spectra. The Lya profile 
can be very complex, consisting of both emission and absorption 
features and this combination often leads to asymmetric profiles 
with a si gnificant amount of absorption blue-w ards of the emis- 
sion line ( Shapl ev et al.l2003l : lKomei et alj201(ll) . For the purposes 
of producing composite spectra of the LBGs, we take a relatively 
simple approach to the measurement of the equivalent widths of 
our galaxy sample. For a given spectrum, we measure an equiva- 
lent width for the emission line if clearly identifiable and if not we 
make a measurement of the absorption profile. To do this, we fit 



The VLT LBG Redshift Survey I 1 3 



Table 5. Redshift ranges of z > 2 galaxies identified from each of our photometric selections. 



Field 



LBOJ^Rll 



LBG.PRI2 



LBGJ'RD LBG.DROP 



Q0042-2627 
JO 124+0044 
HE0940-1050 
J1201+0116 
PKS2126-158 



2.74 ±0.28 



3.02 ±0.33 



2.71 ±0.29 2.45 ±0.41 
2.98 ±0.29 2.72 ±0.27 



2.66 ±0.26 2.67 ±0.30 3.04 ± 0.28 

2.86 ± 0.34 

2.67 ±0.29 2.85 ±0.39 3.10 ±0.21 
2.61 ±0.29 2.74 ±0.33 

n/a 3.30 ±0.29 



All fields 



2.90 ±0.32 2.66 ±0.28 2.67 ±0.30 



I ±0.36 



Q0042-2627 




4.0 -26.52 



Dec, 
-26.19 



4.0 C.47 



4.0 -11-3 



4.0 1-00 



-25.87 

. ^11.4 
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Figure 11. Distribution in R.A., Declination and redshift for each of our five fields. Spectroscopically confirmed LBGs are marked by blue filled circles 
and known QSOs by dark red stars. We also identify those QSOs with low-resolution spectra available (red circles, i.e. VLT VIMOS and AAT AAOmega), 
medium-resolution spectra (red crosses, i.e. SDSS - SDSS J1201-I-01 16 only) and high-resolution spectra (red squares, i.e. VLT UVES, Keck HIRES). 



a polynomial to the continuum and a Gaussian fit to the Lya line 
profile and estimate the equivalent width from these fits. 

The individual LBG spectra were normalized prior to con- 
structing the composite, using the median of the rest-frame UV 
continuum in the range 1300A < Area* < 1500A. After this nor- 
malization, we rescale the LBG spectra to the rest-frame and re- 
binned the spectra before combining the samples to produce the 



final composite spectra. We note that all the spectra were calibrated 
using the VIMOS master response curves prior to this process. 

The composite spectra are shown in Fig. [15] and are split into 
(from bottom to top) equivalent width ranges of W< — 20A (50 
galaxies), -20A<W<0A (134 galaxies), 0A<W<5A (166 galax- 
ies), 5A<W<10A (218 galaxies), 10A<W<20A (181 galaxies), 
20A<W<50A (112 galaxies) and W>50A (60 galaxies). Between 
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Figure 12. Number counts as a function of Rvoga magnitude for all fields, 
except for JO 124+0044 in which Ivcga is used. The shaded histograms show 
the numbers of successfully identified objects with the colour coding the 
same as in Fig. 1101 the cyan histogram shows counts of LBG_DROP objects, 
the red shows LBG_PRI1 objects, the green shows LBG_PRI2 objects and 
the blue shows LBG_PRI1 objects. The unshaded histogram shows the total 
number of candidates observed with VLT-VIMOS in each field (i.e. the gap 
between the shaded regions and solid line shows the number of unidentified 
objects as a function of magnitude). Contamination levels from stars and 
low-redshift galaxies for each field are given by the dashed line in each 
panel. 
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Figure 13. Sky densities of the LBG sample as a function of Rvoga 
magnitude. The red open circles give the total densities of objects in our 
LBG JRIl, LBG JRn and LBG.DROP photometric selections. The 'VLT' 
densities (filled red circles) are estimated using the total photometric den- 
sities multiplied by the fraction of successfully identified LBGs from the 
VLT spectroscopic observations and are corrected for incompleteness in 
the spectroscopic sample at faint magnitudes. Raw (open blue triangles) 
and correcte d (fille d blue triangles) number cou nts are also shown from 
ISteidel et al.l J2003h . Note that we transform the ISteidel et alj J2003l) AB 
system V, magnitudes by —0.14 to convert to -Rvcga ISteidel & Hamilton! 
[i993.1 



I VLT da:a 

I 

I ' Shapley et al. 03 



^ 



1000 2000 

Av(em - abs) (km/s) 



them, the composites incorporate a total of 921 of the galaxy sam- 
ple, excluding any objects with q < 0.5 or with significant con- 
tamination, for example from zeroth order overlap. The key emis- 
sion and absorption features are marked and we can immediately 
identify both absorption and weak emission for the ISM lines: Sill, 
Ol-l-Sill, CII, SilV and CIV. All the features have been marked at 
z — 0. The offset between the line centres of the Lya emission and 
the ISM absorption lines is evident in these composite spectra, a 
result of the asymmetry of the Lya, potentially combined with an 
intrinsic difference between the velocities of the sources of the Lya 
emission and the ISM absorption features. 



Figure 14. Distribution of the velocity offsets between ISM absorption lines 
and the Lya emission line in individual galaxies from our redshift sur- 
vey (solid histogram). We measure a mean velocity offset between Lya 
emission and the ISM lines of AV = 625 ± 510 kms~ . The result of 
IShaplev et alj 120031) , which has a mean of 650 kms~^ is shown by the 
dashed histogram. 



3.6 VLT AGN and QSO observations 

As discussed earlier, we also targeted a small number of z ~ 3 
QSO candidates selected from our UBR photometry. In combina- 
tion with this, due to the similarity in the shape of the spectra of 
LBGs and QSOs, the LBG selections also produced a handful of 
faint QSOs and AGN. We present the spectra of these in Fig. [16] 
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Figure 16. Spectra of the z > 2 QSOs observed as part of the VLT VIMOS LBG survey. Redshifts and ii-band magnitudes are given for each QSO and 
significant broad emission features are marlced. 



whilst the numbers of QSOs in each field are given in tableU The 
positions of the observed QSOs are also shown in Fig. 1111 



4 CLUSTERING 

In this section we present the clustering analysis of the z ~ 
3 galaxy sample, incorporating estimates of the angular auto- 
correlation function for our complete LBG candidates catalogue 
and the redshift space auto-correlation function of our spectro- 
scopically confirmed sample. Developing from these estimates, we 
use a combined sam ple of the VLT VIMOS LBG data-set and the 
ISteidel et alj ( 120031) data-set to evaluate the 2-D correlation func- 
tion and place constraints on the infall parameter, /3, and the bias 
paremeter, b. Finally, we relate the clustering properties of the 
2 ~ 3 sample to those of lower-redshift samples. 



4.1 Angular Auto-correlatlon Function 

We now evaluate the clustering properties of our candidate and 
spectroscopically confirmed LBGs. Using all five of our imaging 



fields, we begin by calculating the angular correlation function of 
the LBG candidates. We use all LBG candidates selected using 
the LBG_PRI1, LBG_PRI2, LBGJRB and LBG.DROP selections 
plus the candidates from the J0124+0044 field. The total number 
of objects is thus 18,489 across an area of l.Sdeg . First we create 
an artificial galaxy catalogue consisting of a randomly generated 
spatial distribution of points within the fields. The angular auto- 
correlation function is then given by the Landy-Szalay estimator 
JLandv & Szaiavlll993l) : 



w{e) = 



(DP) - 2 (PR) + (RR) 
(RR) 



(4) 



where PP is the number of galaxy-galaxy pairs at a given sep- 
aration, 9, PR is the number of galaxy-random pairs and RR is 
the number of random-random pairs. The random catalogues were 
produced within identical fields of view to the data and with sky 
densities of 100 x the real object sky densities, in order to make 
the noise contribution from the random catalogue negligible. We 
estimated the statistical errors on the 1^(9) measurement using the 
jack-knife estimator. 
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Figure 15. Composite spectra collated from our VLT VIMOS sample. Each 
spectrum shows the composite of a sub-sample of the LBGs, grouped by 
LyQ equivalent width measurements. The key UV spectral features dis- 
cussed in the text (i.e. Lya and Ly/3 emission/absorption, ISM absorption 
lines) are all evident in these composite spectra. 



Measurements of w{0) in small fiel ds are subject to a bias 
known as the integral constra int (e.g. iGroth & PeeblesI 119771 : 
|Peebleslll98dlRocheetal.lll993h . This is given by: 



(5) 



(6) 



— / / w{9)dnidn2 



where the 'true' w{9) is then: 



W{9) = (WmeasiO)) + O" 



where {wmeas{d)) is the measured correlation function, averaged 
across the observed fields, and 1^(6) is the correct correlation func- 
tion. As in lRoche et al.l J2002h . we evaluate the integral constraint 
using the numbers of random-random pairs in our fields: 
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Figure 17. The angular correlation function, w{6), from our imaging fields. 
The open stars show the correlation function for the photometrically se- 
lected sample, whilst the filled stars show the same correlation function 
connected for stellar and z < 2 galaxy contamination as described in the 
text. The dashed red and solid red lines show the double power law models 
fitted to the raw and contamination corrected correlation functions respec- 
tively. We also show a m odel determined from the rg, 7 measurements of 
Ida Angela et alj I2005ah - dash-dot blue line. The blue triangles and dot- 
ted line show the correlation function and best fitting power law model for 
the photometrically selected z < 2 galaxy population. The blue dash line 
gives the result of ,Adelberger et al., l ,2Q05b,) . with rp = 4.0 /i~^Mpc and 
7 = 1.57. 



^A 



j:NRR{e)e- 



(7) 



The results of the w{9) calculation for the full photometrically 
selected LBG sample are shown in Fig.[T7](open red stars). 

Additionally we show the correlation function, estimated in 
the same way, for the remaining 23 < i? < 25.5 galaxy population 
(i.e. all galaxies in the given magnitude range not selected by the 
LBG colour selection - blue triangles). This gives an estimate of the 
clustering for the z < 2 galaxy population in the LBG fields. Based 
on the spectroscopic results, we estimate that 60% of the photomet- 
ric selection consists of 2: > 2 galaxies whilst the remaining 40% 
consists of contaminant z < 2 galaxies and galactic stars. In order 
to determine more accurately the clustering of our selected z > 2 
galaxy population, we therefore correct the w{0) measurement for 
the effects of contamination. The correction is given by: 



3(6') = w^<2(6i)/|<2 + ■u;lbg(^)/lbg 



(8) 



where Wmeas is the total measured correlation function, 'Wz<2(0) 
is the correlation function of the contaminant galaxies, fz<2 is the 
fraction of contaminant galaxies, wlbg {d) is the correlation func- 
tion of the z > 2 galaxies and /leg is the fraction of z > 2 
galaxies. We therefore use the measured correlation function (i.e. 
open red stars in Fig.ll7t and the measured z <2 correlation func- 
tion (i.e. blue triangles in Fig.|17|l along with the spectroscopically 
measured fractions of 2 > 2 and z < 2 galaxies to estimate the 
z > 2 galaxy correlation function (i.e. wlbg)- The result is shown 
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by the filled red stars in Fig. [T7] At all scales we find a higher 
measurement of the z > 2 correlation function after applying this 
correction. We note that the wlbg(6') measurement shows signs 
of a change in slope at ~ 0.6 — 1', suggestive of the combina- 
tion of one and two halo terms used in Halo Occupation Distribu- 
tion mode ling (HO D. e.g.lAbazaiian et al.l2005l : IZheng et alj2005l : 
IWake et a l. 2008; Z heng et alj2009l) . 

We now quantify the clustering amplitude of the raw and cor- 
rected w[6) measurements using a simple power-law fit, with con- 
stants Aw and 5 such that: 



w{e) = Ayje- 



(9) 



Fitting to the data to the large scale clustering (0.8' < 9 < 
10') for the uncorrected ■w{6) we obtain best fit parameters of 
Aw = 1.08 ± 0.27 X IQ-^'deg* and S = O.rQtoH- Using the 
same angular range with the corrected w{0) gives parameters of 
Aw = l.85t°Qt\ X 10"^deg* and S = 0.82to:}2- We also perform 
a fit to the z < 2 correlation function. In this case, the cluster- 
ing is fit by a power law with Aw = 2.3lto;gg x lO^^deg* and 
5 = 0.57t°;;ji (dotted blue line in Fig.[T7ll. 

We now estimate the real-space correlation function, 
£(r), from our meas urement of w{8) using Limber's formula 
jPhil lipps et al. 1978) with our measured redshift distribution 
(Fig. not . This is performed for both the raw w{6) and the 
contamination-corrected ■w{6) with a double power-law form of 
^(r) given by: 



,.^{^y\r<r.) 



^2 = 



ro,2 
r 



{r > Tb) 



(10) 



(11) 



where rj is the break at which the power-law is split between 
the two power-laws, ro is the clustering length and 7 is the 
slope (which is given by 7 — 1 + 5). We perform x^ fitting 
over the ro-7 parameter space to both the uncorrected and cor- 
rected w{6) results. Firstly for the uncorrected result, we find 
ro,2 = 3.14t;j:3e /I'^Mpc and 72 = 1.8lto;?4. For the cor- 
rected w{6), we determine a clustering length above the break of 
ro,2 = 4.37l!5:5g /i-^Mpc, with a slope of 72 = 1.61 ± 0.15. 
The full results are given in table|6]and the best-fitting 'w{9) mod- 
els are plotted in Fig.[T7] We note that for continuity in the double 
power-law function, the break is found to be at rj, ~ 1.5 fc~^Mpc. 
Co mparing our result to previous results, Ida Angela et al.l 
( l2005ah obtained a clustering length of rp = 4.481° ?i /t'^ Mpc 
with a slope of 7 = 1.76t|5[5^ and I Adelberger et al.l (l2003h ob- 
tained ro = 3.96 ± 0.15 /i"%pc and 7 = 1.55 ± 0.29, both us- 
ing a single powe r-law function fit (^ (r) = {r/ro)~~') to the same 
2 ~ 3 LBG data JSteidel et alj|2003l) . Our sample appears to have 
a comparable clustering strength, which is slightly higher when 
corrected for stellar/low-redshift galaxy contamination. A further 
comparison can be made with the work of iFoucaud et al.l ( l2003h . 
who measured an amplitude of ro — 5.9 ± 0.5 /i~^Mpc from the 
w{9) of a sample of 1294 20.0 < -Rab < 24.5 LBG candidate s 
in the Canada-France Deep Fields Survey (McCrackenetaLyOOlJ). 
[Hildebrandt et al. (2007) measure the clustering of LBGs in the Ga- 
BoDS data and find a clustering length of ro — 4.8 ± 0.3 /i~^Mpc 
for a s ample of 22.5 < -Rvega < 25.5 galaxies. Subsequently 
to this, iHildebrandt et alj ( 120090 measured the clustering proper- 
ties of LBGs selected in the ugr filters from the CFHTLS data 



and measured a clustering length of ro = 4.25 ± 0.13 /i~^Mpc 
with a magnitude limit of tab < 25 and using redshift estimate s 
based on the HYPERZ photometric code ( Bolzonella et alJuOOOh . 
Our contamination-corrected result appears consis tent with most 

S irevio us work, although lower than the result of iFoucaud et al.l 
2OO3I) . 



4.1.1 Slit Collisions 

After calculating the angular correlation function, we next use the 
redshift information from our spectroscopic survey in order to con- 
firm the clustering properties of the LBGs. However, before we 
do this we need to evaluate the extent to which we are limited in 
observing close-pairs by the VIMOS instrument set up. With the 
LR_Blue grism, each dispersed spectrum covers a length of 570 
pixels on the CCD. Further to this each slit has a length (perpendic- 
ular to the dispersion axis) in the range of 40-120 pixels. Given the 
VIMOS camera pixel scale of 0.205" /pixel, each observed object 
therefore covers a minimum region of ~ 120" x 8.2", in which no 
other object can be targeted. 

In order to evaluate this effect, we calculate the angular auto- 
correlation function for only those candidate objects that were tar- 
geted in our spectroscopic survey, WBiits{d)- To do so we require 
a tailored random catalogue that accounts for the geometry of the 
VIMOS CCD layout. We therefore create random catalogues for 
each sub-field using a mask based on the layout of the four VIMOS 
quadrants, excluding any objects that fall within the 2' gaps be- 
tween adjacent CCDs. The sky-density of randoms in each sub-field 
is set to be 20 x the sky-density of data points in the corresponding 
parent field. From this subset, which consists of ~ 3400 targeted 
objects, we calculate Wsiits(9) using the Landy-Szalay estimator 
(equation |4l(. The ratio of 1 + Wsiit{9) to the original measurement 
of 1 + w{9) (prior to correction for contamination) is shown in 
Fig.[T8](open circles). At S > 2' the two correlation functions fol- 
low each other closely and give a ratio of ~ 1. However at separa- 
tions of ^ < 2' we see an increasingly significant loss of clustering 
showing the effect of the instrument setup. At redshifts of 2 « 3, 
the 2' threshold of the effect corresponds to a comoving separation 
ofr'^ 2.6/i"^Mpc. 

The dashed line in Fig. [TS] shows a fit to the ratio between the 
slit-affected clustering measurement and the original measurement. 
We use this fit to provide a weighting factor dependent on angular 
separation, WsHt{9), which is given by: 



Wsm(9) 



1 



0.07385- 



(12) 



Applying this weighting function to DD pairs at separations of 
9 <2' then allows the recovery of the original correlation function 
from the VIMOS sub-sample correlation function down to separa- 
tions of 5 ~ 0.1'. Below 5 ~ 0.1' however we are unable to recre- 
ate the original candidate correlation function as no close pairs can 
be observed below this scale due to the slit lengths (8" < 9 < 24") 
used in the VIMOS masks. 



4.2 Semi-Projected Correlation Function, Wp (a) 

We next present the semi-projected correlation function Wp{a) for 
the 1020 q ^ 0.5 VLT LBGs. Here, a is the transverse separation 
given by the separation on the sky, whilst -k will be its orthogonal, 
line-of-sight compo nent. We first estimate Wp{a) for the full VLT 
LBG sample using dDavis & Peebleslll983h : 
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Table 6. Clustering results based on the raw w(9) and the w{8) coiTected for stellar and low-redshift galaxy contamination. 



(xlO-^deg"*) 



(/i-lMpc) 



71 



'■0,2 

(/i-^Mpc) 



72 



Uncorrected data 
Contamination-coiTected 



1.08 
1.85 



+0.27 
0.27 

+0.41 
0.21 



0.76 
0.82 



0.07 

0.17 

+0.11 

0.12 



2.16 

3.04 



+0.24 
0.30 

+0.33 
0.34 



-2.49in^ 

-2.48l»;l? 



2.69 
4.37; 



0.20 

0.26 
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-1.61 ±0.15 
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Figure 18. Effect of 'slit collisions' on the measurement of the angular cor- 
relation function, wg . We show the ratio between the clustering of the en- 
tire photometric sample, given by 1 + w{9), and the clustering measured 
from only those objects that have been spectroscopically observed using 
VLTVIMOS, 1 + Wgiits{0). The observational constraints incurred due to 
the constraint of preventing the dispersed spectra from overlapping on the 
instrument CCD lead to a significant reduction in the clustering measure- 
ment at 9 < 2'. The dashed line shows our parameter fit (equation II 2) to 
the measured ratio, which we use to correct subsequent clustering measure- 
ments made using the spectroscopic galaxy sample. 
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Figure 19. Projected correlation function, Wp{a) of the full VLT, Keck 
tSteidel et al...2003.) and the combined samples. The blue dash-triple-dot 
line represents our best-fit (with 7 = 1.8) to the Keck data of rg = 
4.20to Jr 'i~'^Mpc. The dot-dashed line represents the best 7 = 1.8 



fit to the VLT sample with ro = 3.67t°-24 h-^Mpc. The solid line 
represents the best 7 = 1.8 fit to the combined VLT+Keck sample with 



ro 



3.98 



+0.23 



h ^Mpc and the dotted line represents the double- 
p ower-law model fi t ted to the VLT w(8). The dashed line gives the result 
of lAdelberger et alJ <2005M . with m = 4.0 /i^^Mpc and7 = 1.57. 



i 



p(a) = 2/ a'^,7v)dn 



(13) 



We perform the integration over the line of sight range from tt = 
to 100 h~^Mpc. This encompasses much of the bulk of the sig- 
nificant signal in the correlation function and performing the cal- 
culation over a range of reasonable limits showed the result to be 
robust. The VLT 10^(0-) is shown in Fig.[T9]with the best fit cluster- 
ing model determined by a x^ fit to the data shown as a dotted line. 
For the projected correlation function a simple power law form of 
^(r) gives: 



Wp{a)/a 



r(f)r(^) 



r(i) 



(14) 



where r() is the Gamma function. We perform the fit to the data 
using a fixed value for the slope of the function of 7 = 1.8. With 
this value, we obtain ro = 3.67lo:24 h'^Mpc for the full VLT 
sample. Comparing to the initial estimate from the w{9) measure- 
ment in Fig.[T8l we find the Wp(o") measurement gives a somewhat 
lower value for ro. The difference is at the < 2a level and given 



the level of contamination in the photometric sample, we expect the 
Wp{a) measurement to be the more reliable. 

We next compare the VLT result to the LBG Keck sample 
of ISteidel et"aLl ( l2003h . This sample consists of 940 LBGs in the 
redshift range 2.0 < 2 < 3.9, with a mean redshift of {2) — 
2.96 ± 0.29 (compared to 2.0 < 2 < 4.0 and {2) = 2.87 ± 0.34 
for the VLT LBG survey). The survey is based on observations 
within 17 individually observed fields, with most of these being 
~ 8' X 8' with a few exceptions (the largest field being « 15' x 15')- 
The Keck spectroscopic data covers a total area of 0.38deg^, with 
just a small number of the fields being adjacent. The median rest- 
frame UV absolute magnitude is Afi7oo = —17.92 ± 0.02, based 
on the commonly used transformati ons to Mi 700 using the ob- 
served magnitude s 7?. and G (e.g. ISawicki & Thompson! uOOa 
iReddv et al.ll2008l) . With the same method (and the transforma- 
tions to 7?. and G AB magnitudes given by Steidel & Hamilton! 
1 1993!) , we estimate a median rest-frame UV absolute magnitude 
of MnoQ — —18.19 ± 0.03 for our VLT sample. The samples ap- 
pear broadly compatible, with the Keck sample having a marginally 
fainter average absolute magnitude, most likely due to the greater 
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number of fainter objects (7?. > 25) observed with the deeper spec- 
troscopy obtained for the Keck sample. 

Combining the two spectroscopic data-sets gives a total of 
1,980 LBGs over a total area of 1.56deg^. In Fig. [19] we further 
present the Keck and combined results for 'Wp{a). The VLT re- 
sults are slightly lower than for the Keck data in the range 1 < 
a < 7 h~^Mpc. The result for the combined sample is domi- 
nated by Keck pairs for a < 7 h~^Mpc and VLT pairs at larger 
scales. The solid line represents the ro — 4.201q ^5 h~^Mpc, 
7 = 1.8 fit for the Keck data. The dashed line represents ro = 
S-SStaii h'^Mpc, which gives the best 7 = 1.8 fit to the 
VLT+Keck combined Wp data. Also shown is the best 7 = 1.8 
fit to the full VLT sample with ro = 3.67toil h'^Mpc. 

To calculate Wp{a) for the double power-law ^(r) that we fit- 
ted above to the VLT w{6) we used the relation 



Wp{a) 



J a 



r£,(r) 



^/P 



zdr 



(15) 



The dot-dashed line in Fig.[T9]then shows that this model also gives 
a good fit to the combined Wp{a). 



4.3 Redshlft-Space Correlation Function 

The redshift-space correlation function, ^(s), is an estimator of 
the clustering of a galaxy population as a function of the redshift- 
space distance, s, which is given by s = Vo"^ + 1"^- Now, us- 
ing the full VLT sample of 1,020 g ^ 0.5 spectroscopically con- 
firmed z > 2 galaxies, we estimate ^(s) using the simple esti- 
mator ^(s) = DD{s)/DR{s) — 1. Again the random catalogues 
were produced individually for each field to match the VIMOS 
geometry and with 20 x the number of objects as in the associ- 
ated data catalogues. The DD pairs were then corrected for slit 
collisions using the angular weighting function (equation 1 12t ap- 
plied to pairs with separations of ^ < 2'. The result is shown in 
Fig. |20] (filled circles) with Poisson error estimates. The accuracy 
of these errors is supported by analysis of moc k catalogu es gener- 
ated f rom N-body simulations ( Ida Angela et al. 2005a; Ho vle et al.l 
2OO0I). Plotted for com parison is the Keck result as analysed by 
da Angela et al.l OOOSah . Also shown is the combined VLT+Keck 
^(s) result. 

The VLT and Keck samples show good agreement at separa- 
tions of s > 8 h~^Mpc, however the VLT sample shows a signif- 
icant drop in clustering strength at 1 < s < 8 h~^Mpc compared 
to the Keck measurement. This seems at odds with the w{9) result, 
which points to the two samples having similar clustering strengths. 
However, we note that the estimate of the line-of-sight distances is 
sensitive to any intrinsic peculiar velocities and also errors on the 
redshift estimate, which will have a consequent effect on the mea- 
sured redshift space correlation function. In addition to this, the pe- 
culiar velocities are an important element in the cross-correlation 
between the galaxy population and the Lygforest, wh ich is pre- 
sented with this galaxy sample in lCrighton et al.l ( 120101) . We there- 
fore now estimate the effect of our redshift errors on this re- 
sult. The error on a given LBG redshift is a combination of 
the mean error on the spectral feature measurements, which is 
given by the measurement error on the Lya emission line from 
Fig- 13 (i-e- ~ ±450 kms^^ given average spectral S/N=5.5 in 
the full VLT sample) combined with the error on the estimation 
of the redshift from the measurement of the outflow features (~ 
±200 kms"'^). In addition, there will be some contribution from 
intrinsic peculiar velocities. We estimate this contribution based 
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Figure 20. Redshift-space clustering function, ^{s), calculated from 1020 
spectroscopically identified LBGs in the full VLT, Keck and combined 
samples. Also shown is the result from the 529 LBGs in the high S/N 
VLT sample. The models generally adopt the 7 = 1.8, ^(r) ampli- 
tudes fitted to Wp{cr). Thus the combined VLT+Keck model assumes 
ro = 3.98 h^^Mpc and expected velocity dispersions of < w"^ >i/2 = 
720 kms"^ (VLT) and < w^ >i/2= 400 kms'^ (Keck). Also shown 
■ ^ Mpc from the full VLT Wp (a) result and 



is a model with ro 

„2 



3.67 h- 

-1 



< wf >l/2= 1000 kms~^, improving the VLT fit. A further model with 
ro =4.2 /i~^Mpc from the Keck Wp{a) gives a good fit to the Keck g(s) 



with < w)2 >l/2= 600 kms"^ 
w{6) model, assuming < w^ > 
/3 = 0.48 (see Sect.l4!4l) 



Finally, we show the 2-power-law VLT 



1/2= 1000 kms" 



AU models assume 



on the w ork of dTummuangpak et aljin predl . lTummuangpak et al.l 
din preph use the Ga laxies-Intergalactic Medium Interactions Cal- 
culation (GIMIC, ICrain et alj 120090 . which sa mples a num- 
ber o f sub-grids of the Millennium Simulation ISpringel et al.l 
(|2005|), p opulating these with baryons using hydrodynamic sim- 
ulations. iTummuangpak et al.l ([In prep) measure a mean intrin- 
sic peculiar velocity based on galaxies in the GIMIC simu- 
lations in redshift slices at z = 3.06 and find a value of 
~ 140 kins^^. Combining this in quadrature with the esti- 
mated measurement errors gives an overall velocity dispersion 
of a, = y/(450 kms-^)2 + (200 kms-^)2 + (140 kms~^)2 ^ 
510 kms^^. The expected overall VLT pairwise velocity disper- 
sion is therefore < wl >'^^^= ^ x 510 ^ 720 kms"\ Sub- 
stituting a Lya emission-line velocity error o f ±150 kins~^ ( based 
on a measurement error of A2; ~ 0.002 from ISteidel et al.l20 03) in 
the above expression similarly implies an expected < mf >^"Ri 
400 kms~^ for the Keck pairwise velocity dispersion. 

On small scales, the above random pair-wise velocity disper- 
sion leads to the well known 'finger-of-god' effect on redshift- 
space maps and correlation functions. On larger scales, bulk infall 
motion towards over-dense regions becomes a significant factor and 
causes a flattening in the line-of-sight direction in redshift space. 
We now model these two effects to see if the ^(r) estimates mea- 
sured from the LBG semi-projected correlation function, Wp{a), 
and the angular correlation function, w{9), are consistent with the 
me asured LBG redshift-s pace correlation function, £,{s). Follow- 
ing iHawkins et alj ( |2003|) . we use t he real-space pre scription for 
the large scale infall effects given bv lHamiltonI ( 119921) whereby the 
2-D infall affected correlation function is given by: 
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C'((7,^) = Co (s)Po(m)+ 6(5)^2 (M)+e4(s)P4(M) (16) 

where Pi{^) are Legendre polynomials, ^ — cos(6) and 6 is the 
angle between r and tt. For a simple power-law form of ^(r) the 
forms of Ci(s) are: 



Us) 



1 + 



2/3 , /32 



i{r) 



Us) 



3 



+ 



4^2 






7(2 + 7) 
(3-7)(5-7) 



U) 



ar) 



(17) 



(18) 



(19) 



where 7 is the slope of the power-law form of the real-space cor- 
relation function: ^{r) = (r/ro) '. For the 2-p ower-law model 
case we use th e equivalent ex pressions derived by Ida Angela et alj 
( l2005ah . As in lHawkins et alJ ( [200a) . the infall affected clustering, 
^'{cr, tt) is then convolved with the random motion (in this case the 
pair-wise motion combined with the measurement uncertainties): 



C(<^,^) 



^'{a,n-v{l + z)/H(z))f{v)dv 



(20) 



where H{z) is Hubble's constant at a given redshift, z, and f{v) is 
the profile of the random velocities, v, for which we use a Gaussian 
with width equal to the pair-wise velocity dispersion, < w^ >^ ". 

With this form of f{v), we take the expected pair-wise veloc- 
ity dispersion, < w^ >i/2— 720 kms~^ for the full VLT sample 
and < w^ >i"= 400 kms~^ for the Keck sample. Now taking 
an estimate of /3 = 0.48 (see Section l4!4t . we may model the ef- 
fect of these velocity components on the LBG sample C(o", vr), first 
using the single power-law fit to the combined sample Wp(a) with 
ro — 3.98 h~^Mpc and 7 — 1.8. The form of C(s) estimated 
from the resultant C(o", ti") is plotted in Fig. [20] (solid line). While 
the model with < w^ >i/2= 4qq kms~^ gives a good fit to the 
Keck data, the model with < w^ >i/2_ 720 kms~^ appears to 
overestimate the VLT correlation function at s < 8 h^^Mpc. Even 
increasing the velocity dispersion to 1000 kms~^ did not signifi- 
cantly improve the fit. We also analysed the LBG sub-sample de- 
fined by having spectral S/N > 5. We found that C(s) for this 
subsample did rise and would require a pair-wise velocity disper- 
sion of ~ 1000 kms^^ for the model to fit the data. This is sig- 
nificantly more than the predicted pair-wise velocity dispersion 
of « 600 kms~^, calculated by replacing the velocity error of 
450 kms~^ for the full sample by 350 kms^^ in this case, cor- 
responding to average S/N=8.25 in Fig. [7] The fact that the points 
at s < 1 ^~^Mpc and those at s > 8 h~^Mpc agree with the 
model argues against an even larger velocity dispersion. 

The other possibility is that the ro = 3.98 ft^^Mpc model 
may be too high for the VLT (,{r). Certainly the amplitude of ^{r) 
from the VLT Wp{a) appears lower than either that from the VLT 
w{9) or the Keck Wp{a). Fig. |20] shows that the fit improves for 
the full VLT samples and the high S/N subsample if the correlation 
function amplitude reduces to ro — 3.67 h^^Mpc as fitted to the 
VLT uip(cr), coupled with the velocity dispersion increasing to < 
wl >i/2= 1000kms~\ 

The combined VLT+Keck sample is very similar to the Keck 
sample at small scales. Even for the Keck sample we find that an in- 
creased pairwise velocity dispersion of < w^ >^''^^ 600 kms"^ 



4.2 h- 

1 



^Mpc. Fo r the Keck LBGs, 



JSteidelet all 12 003) + intrinsic 
. lAdelberger et alj|2003f) combines in 



is needed to fit C(s) if ro 
the velocity error (±150 kms 
outflow error (±200 kms 
quadrature to give ±250 kms"^ as the error for the line mea- 
surement. Subtracting from ±600/\/2 kms~^ would imply ~ 
340 kms"^ for the pairwise intrinsic velocity dispersion. Clearly 
for the VLT samples the implied velocity dispersion would be even 
larger. 

We have also used the double power-law £^{r) indicated by the 
VLT w{9) to predict ^{s). Since the steepening takes place at r < 
3 h~^Mpc, this means that we would need even higher velocity 
dispersions to fit S,{s). Fig. 1201 shows that the double power-law 
model needs at least a velocity dispersion of ~ 1000 kms~^ to fit 
the VLT±Keck combined sample. 

We conclude that the low ^{s) we find in the full VLT sam- 
ple may be caused by a statistical fluctuation in the LBG clustering 
due to a lower than average ro and a higher than average veloc- 
ity dispersion. The VLT sample is designed to improve correlation 
function accuracy at large scales, particularly in the angular direc- 
tion, and the somewhat noisy result for C(s) at the smallest scales 
reflects this. Overall, we conclude that the velocity dispersions re- 
quired by C(s) ar e bigger than reported pr eviously for the Keck data 
(400 kms"^ by Ida Angela et"aLll2005al) with the Keck and VLT 



samples now being fitted by < w^ > 



l/2_ 



600 - 1000 kms" 



close to what is expected from estimates of the redshift errors. 



4.4 Estimating the LBG infall parameter, I3{z — 3) 

The infall parameter, /3, quantifies the extent of large scale coherent 
infall towards overdense regions via the imprint of the infall motion 
on the observed redshift space distortions. Given its dependence on 
the distribution of matter, measuring /? can provi de a useful dynam- 
ical constraint on O^ ( z) ^ HamiltonI 19 92: Heave ns & TavloJ 19951: 



iHaw kins et"al] l2003l: Ida Angela eTaLJ 120081 : iCabre & Ga ztaiiagal 
2009). It relates the real-space clustering and redshift-space clus- 
tering as outlined in the previous section (see equations 11 6|to[T9t. 

We shall measure I3{z — 3), usin g the combin a tion o f our 
VLT LBG data and the LBG data of Isteidel et"ai] ( |2003|) . As 
noted above, the VLT and Keck samples complement each other 
in the wide range of separation, a, in the angular direction for 
the VLT sample and the high sky densities of the Keck samples, 
which help define the clustering better at small scales. As dis- 
cussed in section 14.31 the two samples possess comparable real- 
space clustering strengths, with measured clustering lengths of 
ro = 3.67tl^;24 h-^Mpc and ro = 4.20l^:ig h~^Mpc for the 
VLT and Keck LBG samples respectively. The higher estimated 
velocity error of the VLT sample at ±450 kms~^ compared to 
the Keck ±300 kms^^ will make little difference due to the fur- 
ther contributions of the outflow errors and intrinsic velocity dis- 
persions, the dominance of the Keck data at small scales and the 
smaller effect of velocity errors at large spatial scales where the 
VLT data is dominant. We shall therefore combine the two samples 
in the two methods we use to measure /3. 

We first estimate /3 by simply comparing the amplitude of C(s) 
and ^{r) and using equation [TT] at large scales. Fig. 121! shows the 
C(s) from the combined VLT and Keck samples divided by the 
best fit model for C(r) from the semi-projected correlation function, 
Wp{a), with ro = 3.98lo;i2 h'^Mpc and 7 = 1.8. Equation [TTl 
applies only in the linear regime, so we do not expect it to fit at 
small separations. We therefore fit at s > 10 h~^Mpc. Fitting in 
the ranges 10 < s < 25 h'^Mpc and 10 < s < 60 h~^Mpc 
gives the two dashed lines in Fig.|2T] which correspond to I3{z — 
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Figure 21. The redshift space correlation function, g(s) divided by the 
real space correlation function, ^(r), with the latter assuming ro = 
3.98 /i~^Mpc and 7 = 1.8. The short and long dashed Unes repre- 
sent the best fit to the data in the ranges 10 < s < 25 /i~^Mpc and 
10 < s < 60 h-iMpc, which correspond to I3{z = 3) = 0.51^023 and 
fS{z = 3) = 0.381 



"^ -0 23 from eQuation fm 



0.38t|5;23 with the difference 



3) = 0.51+^:23 and l3{z = 3) 

between these two giving a further estimate of the uncertainty in /3 

from this method. 

We next estimate /3 using the shape of the 2-point correlation 
function, ^(a, n), to measure the effect of redshift space distortions. 
We calculate ^(a, n) for the combined sample. As with our de- 
termination of ^(s), we use the simple DD/DR estimator taking 
randoms tailored to each individual field, with errors again calcu- 
lated using the Poisson estimate. The resultant ^{a, tt) is plotted in 
Fig. |22] The elongation in the n dimension, due to the pair-wise 
velocity dispersion and redshift errors, is clearly evident at small 
scales. 

Now using this measurement of ^{a, tt), we make an estimate 
of the infall parameter, /?. For this we use the single power-law 
model of ^{r) with ro = 3.98 h'^Mpc and 7 — 1.8 based on 
the semi-projected correlation function of the combined data in 
Fig. [19] With these parameters set, we calculate the model outlined 
in equations 1161 to l2Ql over a range of values of < wf y^''^ and 
/3. We then perform a simple Ax^ fitting analysis and jointly es- 



timate < w^ > 



l/2_ 



700 ± 100 kms~ and infall parameter of 
/3lbg(-z « 3) = 0.48 ± 0.17 for our combined LBG sample. The 
contour plot of Ax^ for the fit in the < tu^ >^ ": /3 plane is given 
in Fig.[S] 

We note that if we allow the amplitude of ^(r) to be fitted 
as well as the other two parameters, then the results move to /3 = 
l.l±0.4and< w^ >'^^^= 800±100 kms"^ forabestfit7 = 1.8 
value of ro — 3.64 h~^Mpc. Taking the Keck sample on its own, 
we again find /? = 0.9 - 1.5 and < wl >i/2= 650 - 750 kms^^ 
if ro is not or is allowed to float respectively. The Keck fits have 
to be resticted to s < 25 h^^Mpc because of the small a range in 
the angular direction and if we apply the same cut to the combined 
sample, values of /3 again rise to /? = 0.8 — 1.1 and < w^ >^"Ri 
800 kms^^, similar to the results for the Keck sample. Although 



Figure 22. ^(ct, tt) projected correlation function calculated from the s pec- 
troscopically confirmed LBGs from the combined ISteidel et alj 42003h and 
VLT VIMOS LBG samples. The best fit model contours are marked as solid 
lines with /3(2 = 3) = 0.48 and < w^ >l/2= 700 kms-^ 



the errors are clearly still significant, we prefer values of /? ~ 0.5 — 
0.6 given by the amplitude of ^(s) and the shape of ^(a, tt) for the 
combined sample which seems best to exploit the advantages of the 
Keck sample at small scales and the VLT sample at large scales. 

We have also checked the effect of assuming the double 
power-law model fitted to the LBG w{0) in Fig. [17] with ro.i — 
3.19 h-^Mpc, 71 = 2.45, ro,2 = 4.37 ft-^Mpc, 72 = 1.61 
and Tb = 1 h~^Mpc. The best ^(cr, tt) fits are then given by 



P 



0.20 ± 0.2 and < 



> 



^^^= 750 ± 150 kms"^ The 



reduced x^ was 3.44 compared to 3.16 for the single power-law 
model. However, allowing the ^(r) amplitude to vary gave /3 — 
0.48lJ!;33 and < wl >i/2= 725tlll kms"^ with fitted ampfi- 
tudes ~ 80% below those estimated from w{9). The small scale 
rise at r < 1 h^^Mpc will not affect our fit much because of the 
lack of statistical power at small separations. Also, the models we 
are using are expected to be accurate only in the linear regime at 
larger scales. The 80% reduction of the amplitude to the large scale 
power-law implies an ro — 4.05 h^^Mpc which is close to the 
ro = 3.98 h^^Mpc value assumed for our single power-law fits 
above, leading to similar fitted values for /3 and < w^ >i/2 ju 
these two cases. The lower /3 from the actual 2 power-law model is 
simply a result of the high ^(r) amplitude implied by w{0) forcing 
P down in the ^{a, n) fit according to equationfTTl 

Comparing our result of 13 — 0.48 ± 0.17 to previous esti- 
mates of I3{z ~ 3), we generally find somewhat higher values than 
da Angela et al. (2005a), who estimate a value of /? = 0.15lo:?5- 
This is partly because we have assumed Qm {z = 0) = 0.3 and fit- 
ted for the velocity dispersion < w'^ >i/2 whereas lda Angela et al.l 
l l2005ij) assumed < w^ >'^/^= 400 kms"^ and fitted for fJ^ (z = 
0). If we assume < w^ >i/2^ 499 kms"^ for the VLT + 
Keck samples, our estimate of /3 reduces to /? = 0.18 for the 
combined sample. The assumption of < w'^ >i/2_ 499 kjiig^i 
seems to be the main factor that drove /3 to lo wer values, also 
helped by the different model for ^(r) assumed bv lda Angela et al.l 
( l2005ah . a 2-power law model with 71 — 1.3 and 72 = 3.29 with 
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Figure 23. LBG pairwise velocity dispersion (< w^ >^/^)-infall 
parameter(/3) Ax^ contours for the VLT+Keck sample, fitting to ^{u, tt) 
< 40 h-'^Mpc. The best fit values are [3 = 0.48 ± 0.17 and 



with 
< w] 

7 = '- 



>l/2 = 



700 ± 100 kms 



-1 



assuming ro 



3.98 h-'^Mpc and 



^6 



— 9 h '^Mpc motivated by fitting the form of ^(s). The con- 



tours in the < w^ > ' : P plane in Fig. |23] show that /3 and 
< w'j. >^" are degenerate - higher j3 implies more flattening 
in the tt direction which can be counteracted by fitting a higher 



< w^ > 



1/2 



to produce elongation in tt. A flatter small scale slope 
for 5(r) also allows a smaller < w1 >^" to be fitted which can 
then allow lower values of /? to be fitted. We have also fitted our 
combined data with a further 2-power-law form for ^(r), now with 
ro,i = 3.98 h-^Mpc, 71 = 1.8, ro,2 = 5.99 h'^Mpc, 72 = 2.6 
and rj, = 15 /i~^Mpc but we find that the results for < wl >^" 
and P from the combined sample are similar to those for the single 
power-law model. 

As well as the higher value of /?, we note that we are also 
fitting higher velocity dispersion values to the combined sample. 



Again the degeneracy between < 



> 



1/2 



and P may be the 



cause. However, the need for high velocity dispersions was also 
noted in the small scale fits to ^(s) particularly for the VLT sample 
but also for the Keck sample. Even < wl >^"= 600 kms^^ 
for the Keck sample implies an intrinsic velocity dispersion of 
< wl >^' ^^ 440 kms^^ taking into account velocity and outflow 



errors on the redshift, much higher than < w^ > 



l/2_ 



200 kms" 



expected from the simulations. If our velocity errors were underes- 
timated then this could be a cause but they would have to be under- 
estimated in both the Keck and VLT datasets. Larger velocity er- 
rors are also contradicted by the consistent widths of the emission- 
absorption difference histograms in Fig. [14] For example, assum- 
ing ±450 kms^^ for the VLT emission velocity error is consistent 
with ±200 kms~^ for the outflow error and ±130 kms"^ for the 
absorption line error. 

We conclude that for ^m{z = 0) = 0.3, the combined 
survey is best fitted by < wl >'^/^= 700 ± 100 kms~^ with 
P = 0.48 ± 0.17 for a single power-law model with 7 = 1.8 
and ro = 3.98 /i"^Mpc. Based on the /3 = 0.49 ± 0.09 value, 
ro = 5.05 /t~^Mp c and 7 = 1.8 values found for 2dFGRS 
JHawkins et alj|2003h linear theory predicts P{z = 3) = 0.22 in 
the rJm = 1 case and P = 0.37 in the Q.m (z = 0) = 0.3 case, with 



ro — 3.98 /i~^Mpc for the latter and transformed appropriately for 
Q.m = 1. Our measurements appear to produce values of P that are 
marginally more acceptable with Q,m{Q) = 0.3 than fim(O) = 1 
but neither case is rejected at high significance; P = 0.22 with 
< wl >^/2— 500 kms^^ is rejected only at 1.5(7 in Fig. [23] More 
importantly, these measurements provide a useful check of the im- 
pact of small- and large-scale dynamics on our measurement of the 
clustering of our 2 ~ 3 galaxies. The estimates of < w'i >^' •^ will 
also be useful in interpreting the effect of star-formation feedback 
from our LBGs on the IGM as measured by the Lyman-alpha forest 
in background QSOs JCrighton et alj2010l) . 



4.5 Estimating the LBG bias parameter, b{z = 3) 

We can now estimate the bias, b, of the VLT+Keck LBG sample 
from our P measurements. The bias gives the relationship between 
the galaxy clustering and the underlying dark matter clustering: 



^n 



b% 



(21) 



where ^j^ij^j is the volume averaged clustering of the dark matter 
distribution and ^ is the volume averaged clustering of a given 
galaxy distribution. In a spatially flat universe, the relationship be- 
tw een the bias, b, and the infall parameter, /3, can be approximated 
bv JLahavetalJl99lh : 



/3 






(22) 



Using this relation with our estimate of P = 0.48 ± 0.17 and 
assuming that Q.m{z = 0) = 0.3 and then given that Q.m{z — 
3) = 0.98, this implies b{z = 3) = 2.06J;o.53- 

We now compare this to an estimate of the bias from our ear- 
lier clustering analysis using equation [21] To do this we calculate 
the dark matter clustering using the CA MB software incorp orating 
the HALOFIT model of non-Unearities JSmith et al]|2003h . From 
this we determine a second estimate of the bias using equation I2l1 
and c alculating the volume averaged clustering function (Peebleg 
1198(1) within a radius, x, for our galaxy sample and the dark mat- 
ter: 



C(^) 



rx 

7 ' 

Jo 



^{r)dr 



(23) 



where ^(r) is the 2-point clustering function as a function of sep- 
aration, r. We use an integration limit of a; = 20 h~^Mpc, en- 
suring a significant signal, whilst still being dominated by linear 
scales. Taking the volume averaged non-linear matter clustering, 
with the volume averaged clustering of our galaxy sample (with 
ro = 3.98 h-^Mpc and 7 = 1.8, from the VLT±Keck Wp{a) 
measurement) and determining the bias using equation [2T] we find 
b — 2.22 ± 0.16, consistent with the estimate from the bulk flow 
measurement of /J = 0.48 ± 0.17 which implies b = 2.06+^53. 
Both values are somewhat lower than the measurement of the bias 
of a sample of LBG s from the Canada-France Deep Survey by 
JFoucaud et alj ( l2003b who measured a value of 6 = 3.5 ± 0.3. 

We now estimate the m ass of typical host haloes for the 
2 ~ 3 LBG sample using the ISheth et al. (2001) prescription for 
the relation between halo mass and bias, determining a host halo 
mass of Mp M = xlO"-^='=°-^/i-~^MQ. Comparing this to o ther 
LBG samples , JFoucaud et alj l l2003l) , lHildebrandt et alj ^2mh and 
JYoshida et al.l yOOSJ) measure halo masses of bright z ~ 3 LBG 
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Figure 24. The shaded regions are the Qf^-j3{z) contours for the 
VLT+Keck sample, fitting to ^{a, tt) with s < 40 ft-^Mpc. The dashed 
lines are the 1- and l-a contours from comparing the z fa 3 LBG and 
the 2dFGRS z ft! 0.1 clustering amplitudes and also using the 2dFGRS 
/3(z ft! 0.1) = 0.49 ± 0.09 result. The dotted Hne is the 1 - cr joint 
contour from applying both of these constraints. 



samples of Mdm ~ lO^^/i~^M0. This difference in mass esti- 
mates reflects the deeper magnitude limits of our survey compared 
to a number of the above results and also a slightly lower red- 
shift range that contribute t o our LBG selection samp ling a lower 
mass range. Work using the Steidel et al.l <200 3112004 ) data, which 
is closer to our own in red shift and depth, report halo masses of 
Mdn. 



lO"-^/i"^M0 JAdelberger etal] l2005bl ; Fconrov et all 
l2008 i). which is closer to the estimate presented here, although our 
result is still somewhat low. 



4.6 Further test of the standard cosmology 

Following the analysis of Ida Angela et al J ( l2005ah we can make 
a further test of the standard cosmology by directly comparing 
the independent values of the bias from the z-space distortion and 
the LBG clustering amplitude. Whereas in the above case we as- 
sumed the DM clustering for the standard model, here we simply 
assume the 2dFGRS clustering scale length which we approximate 
as ro = 5.0 /i~^Mpc and 7 = 1.8. We also assume their value 
of /?(z « 0.1) = 0.49 =b 0.09 from red shift space distortions. In 
similar fashion to lda Angela et al J ( l2005ah we can then for any SlJ^, 
find the mass clustering amplitude at « = 3 and then we can find the 
LBG bias, b{z = 3), by comparing this to the amplitude of LBG 
clustering given by ro — 3.98 /i~^Mpc and 7 = 1.8. This can 
then be converted to /3(z = 3) by using the value for Q,m{z = 3) 
implied by the assumed n°„ and therefore the I3{z) : SlJJ^ relation 
can be drawn. The 1- and 2-a upper and lower limits on this rela- 
tion are shown in Fig.|24] These are overlaid on the Ax^ contours 
(greyscale) from a similar redshift-space distortion analysis as seen 
in Fig.|23]but now allowing Q}^^ and /?(2 = 3) to vary while keep- 
ing < Wj >^'2= 700 kms^^ constant. In this case we have also 
allowed the LBG clustering amplitude to be fitted within a 50% 
range; this is to ensure that the dynamical constraint is as inde- 
pendent as possible of the other constraint which is directly taken 



from the LBG clustering amplitude. We see that although the best 
fit from redshift-space distortions has now moved to lower ft^ and 
lower I3{z = 3), there is still a good overlap between the ±1 — a 
regions of both constraints. The 1 — a joint contours from both 
constraints are shown by the dotted line with the best joint-fit be- 
ing Q,^ = 0.2 and I3{z = 3) = 0.45. Thus there is certainly no 
inconsistency with the standard ACDM model although, as before, 
the 0,%^ — 1 model is still rejected at less than the 2a level. With 
the values of fiJJ^ in a reasonable range, there appears no inconsis- 
tency with the evolution of gravitational growth rates as predicted 
by Einstein gravity, extending the results presented bv lGuzzo et al.l 
fcoOS.) to 2 fti 3. 



4.7 Clustering Evolution 

The space density and clustering evolution of LBGs have fre- 
quently been used to infer their descendant galaxy populations 
at the present day. Initially, their relatively high clustering ampli- 
tudes were taken to mean that they would evolve on standard halo 
models into luminous red galaxies in the richest galaxy clu sters at 
z = (Steidel et al. 1996; Gover nato et al 1998; Adelberg er et al.l 
l2005bh . On the other hand, Met calfe et alj ( 1 19961 , [2001 ) noted that 
the comoving density of LBGs was close to that of local spirals. 
Indeed, they showed that a simple, Bruzual & Chariot (1996), pure 
luminosity evolution model with e-folding time, r = 9Gyr, plus a 
small amount of dust, could explain the LBG luminosity function 
at z ~ 3. Recently, more detailed merger tree models have been 
used to interpret LB G space densities and clustering. For example, 
IConrov et al.l 1200a) have concluded on this basis that the descen- 
dants are varied, with LBGs evolving to become both blue and red 
L* and sub-L* galaxies. 

We now qualitatively compare the clustering strength of our 
LBG samples to that of lower redshift galaxies. We first determine 
the volume-averaged correlation function at 20 h^^Mpc using the 
single power-law form of the clustering of both our own and the 
Keck LBG sample as prescribed in equation [23] The ^(20) mea- 
sured for the VLT LBG sample is shown in Fig. [25] compared 
to a number of measures of the clustering of other galaxy sam- 
ples across a range of redshifts. The VLT+Keck result (ro = 
3.98 /i"^Mpc, 7 = 1.8, z = 2.87) is shown by the filled 
star. We also sho w the measure for the Keck LBG sample alone 
(open star) and the lFoucaud et al.l J2003r) LBG sample (cross). The 
apparent B-band magnitude range of the VLT+Keck sample is 
B = 25.69 + 0.76. Using the overall redshift range of the sam- 
ple (z = 2.87 ± . 34) an d K+e corrections determined using the 
iBruzual & Charlod 120031) stellar population evolution, this equates 
to an absolute B-band magnitude of Mb ~— 21.5 + 1.1. 

For comparison with our data, we have also plotted the esti- 
mated volume-averaged correlation function values for a number 
of low and high redshift galaxy s amples. The open and fi lled red 
triangles show the LRG samples oflSawangwit et al J 120091 ) . giving 
the clustering for a 2L* and 31/* sample respectively (and hav- 
ing absolute i-band magnitudes of A^hab) = —22.4 ± 0.5 and 
Mi(AB) = —22.6 ± 0.4). The open squares show the clustering 
of late -type galaxies from the 2dFGRS as given bv iNorberg et al.l 
120021) with the individual points giving the clustering of galax- 
ies in the absolute magnitude ranges of —18 > Altj > —19, 
-19 > Mtj > -20, -20 > Mbj > -21 and -20.5 > 
Mtj > —21.5 (in order of lowest to highest cl ustering data-points) . 



ita-po 

In addition we plot the blue spiral galaxie s of iBielbv et al.l 20101) 
with the open upside-down triangles and IBlake et alj 120101) with 
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fil led upside-down trian gles, plus the sBzKs (open blue diamond) 
of lHavashietalJ ( l2007h . 

As an illustration of how we may expect the clustering of the 
samples to evolve with time, we first consider a model based on 
the simulated merger history of dark matter haloes (dot-dot-dot- 
dash U ne) calculated from the simulations of iGonzalez & Padillal 
(|2010|), w hilst the method us ed to follow the merger trees is de- 
scribed in lPadillaet alj (12010'). The simulation was performed us- 
ing parameter values of Om = 0.26, SIa = 0.74, ag, — 0.80 and 
Us = 1.0 and consisted of a box size of Lt,ox = 123 h~^Mpc 
containing 512^ particles with a particle mass of lO®/i~^M0. The 
normalization to the LBG data was performed by finding the halo 
mass (lO"-^^*°°*/i~^M0) for which the halo I20 matches the I20 
measurement for the VLT LBG sample at its mean redshift. We see 
that the model predicts little change in the clustering amplitude at 
z = 1 and then stronger evolution to a higher clustering amplitude 
at 2; = 0. The amplitude of the clustering at z=0 app ears consistent 
with that of late- type galaxies in the 2dFGRS survey JNorberg et al.l 
1 20021) . The predicted descendant number density at z = based on 
the halo merger tree model is logjo (n/(/i^Mpc"^)) = -3.491°-^'? 
and is also consistent with number density of the iNorberg et al.l 
( l2002h -20.5 > Mtj > -21.5 late-type population, which is 
equal to logi(,(n/(/!,^Mpc"^)) = -3.64t°;°J. These models are 
able to estimate the transition scale between the 1-halo and 2-halo 
terms in the correlation function of O.TllJgj h^^Mpc, consis- 
tent with the transition scale of rt, ~ 1.5 ± 0.3 h~^Mpc in our 
measured LB G w(^). Overall, the se conclus ions are not dissimi - 
lar to those of IConrov etal] ( |2008|) . However. IConrov et alj J200i) 
predicted higher clustering amplitudes, ro ~ 5 — 6h~^Mpc or 
1(20) = 0.21 - 0.29, at 2 ?5i 1 and ro ~ 6 - 7h"^Mpc or 
1(20) = 0.29 - 0.38, at 2 « for the LBG desce ndants. Given 
these differ ences between the m erger-tree models of IPadilla et al.l 
( 120100 and IConrov et al.l (|2008), we conclude the results appear 
somewhat model dependent. 

We next compare the ^(20) results to simpler clustering mod- 
els. This approach i s partl y motivated by the interpretation of 
iMetcalfe et alj ( 119961 120011) whose passive luminosity evolution 
(PLE) models connected the LBG population at 2 ~ 3 to the late- 
type population at 2; ~ 0. Such models assume that the comoving 
density of the LBG/late-types remains constant with time and the 
clustering models considered here also make this assumption. Al- 
though the models do not take into account halo mergers, it has 
been shown that in the case of Luminous Red Galaxies, such mod- 
els can still provide useful phe nomenological fits to LRG cluster- 
ing out to significant redshifts (^Wake et alJ l2008l : ISawangwit et al.l 
[2009). Therefore we first plot in Fig.l25lthree simple clustering evo- 
lution models: the long-lived model (dashed blue lines), stable clus- 
tering (dot-dashed cyan lines), and no evolution of the comoving- 
space clustering (short-dashed line). All the models have been nor- 
malised to the VLT LBG clustering amplitudes. 

The long-lived model is equivalent to assuming that the galax- 
ies have ages of order the Hubble time. The clustering evolution 
is then governed by their motion within the gravitationa l potential 
and assuming no merging jFrvll 19961 : ICroom et al.ll2005l) . The bias 
evolution is thus governed by: 



6(2) 



1 + 



m - 1 

D{z) 



(24) 



the bias evolution in conjunction with the dark matter clustering 
evolution, again determined using the CAMB software incorporat- 
ing the HALOFIT model. This is then normalized to the measured 
LBG clustering at the appropriate redshift. 

The stable clustering model repre sents the evolut ion of virial- 
ized structures and is characterised by l lPeacockll 19991) : 



ar,^) 



{1 + zy 



(25) 



where r is the comoving distance. 

Finally, the no-evolution model simply assumes that there 
is no evolution of the clustering in comoving coordinates. From 
Eq [24] this model can be thought of as a long-lived model in the 
limit of very high bias, (6(0) >> 1) since then 6(2) ~ b{0)/D{z). 

Evaluating the clustering evolution of the LBGs, first using 
the stable clustering prescription, we would expect the clustering 
of the 2 ~ 3 galaxies to evolve t o a level comparable to that of 
low-redshift LRG galaxy samples JSawangwit et alj|2009l) . giving 
a highly clustered m odem day population. However, as argued by 
IConrov et al.l ( l2008h . the number density of luminous, early-type 
galaxies may not match that of LBGs at 2 ~ 3 as required by this 
virialised clustering model. Alternatively, on the basis of the long- 
lived model, the LBG descendants could either be lower luminos- 
ity red galaxies or higher luminosity blue galaxies. The space den- 
sity of such galaxies is probably more consistent with that of the 
LBG population. This assumes the ACDM cosmology and its spe- 
cific value of as = 0.80. For a lower mass clustering amplitude 
the long-lived model would have higher bias and the 2 = pre- 
dicted amplitude would reduce to more resemble the no-evolution 
model. In this case, the descendants of high redshift LBGs could 
even be the rela tively poorly clustered, star-forming galaxies of 
iBlake et alJ ( l201Q;l . Thus the long-lived models tend to make LBGs 
the progenitors of bluer, or lower-luminosity red galaxies at the 
present day, similar to the conclusion from the merger-tree model 
of Conrov et al. ( 2008). The no-evolution (or long-lived, high bias) 
model would suggest LBGs are the progenitors of bluer galaxies 
with lower clustering ampli tudes, more similar to the conclusions 



of the merger-tree models of lPadilla et alJj201C ) or the simp le pure 



ot the merger-tree models ot IPadiila et alJlzOll ) or the simp 
luminosity evolution models of IMetcalfe et al.l (l99g,|200l|). 



where D{z) is the linear growth rate and is determined using the 
fitting formulae of ICarroll et al.l ( Il992n . We evaluate ^(20) using 



5 CONCLUSIONS 

In this paper we have described the VLT VIMOS survey of 2 « 3 
galaxies in a number of fields around bright 2 > 3 QSOs. In total 
this survey has so far produced a total of 1020 LBGs at redshifts of 
2 < 2 < 3.5 over a total area of l.lSdeg^. This concludes the data 
acquisition for the initial phase of the VLT VIMOS LBG Survey. 
At the time of writing, these are the most up to date observations, 
however the survey has a number of observations only recently ac- 
quired, comprising another 25 VIMOS pointings. Upon comple- 
tion, the survey will comprise a total of 45 VIMOS pointings, build- 
ing significantly on this initial data-set and providing a catalogue of 
~ 2, 000 2 > 2 galaxies over a sky area of 2.11deg^. The wide an- 
gular coverage of VLT VIMOS makes the new LBG study very 
complementary to the previous Keck study which has higher space 
densities over smaller areas and hence increased power at the small- 
est LBG separations but little information in the angular direction 
beyond 10 h~^Mpc. We therefore have frequently used the two 
surveys in combination in the studies of LBG clustering we have 
presented here. 

Based on the fraction of objects observed for this initial VLT 
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Figure 25. The volume-averaged correlation function, ^(20), is plotted 
for our LBG sample alongside 5(20) measurements for several other 
galaxy populations, including L RGs at 2: < 1 (Sa wangwlt et al. 2009), 
star-forming galaxies at 2 < 1 jNorberg et al.ll200Z : .Bielbj^_et_aljj2010|: 
Blake et al. 2010 ) and z ~ 2 and other LBG populations jAdelberger et alj 
2Q03i : iFoucaud et alj2003l) . Further to the observational data, the solid line 
shows the estimated evolu tion of the underly ing dark-matter clustering us- 
ing the CAMB software ('L ewis et al] |2000'). whilst the horizontal dotted 
line, dashed line and dot-dash line show the clustering evolution given no 
evolution in comoving coordinates, the long-lived model and the stable 
model. The dot-dot-dot-dash line shows the clustering evolution based on 
the modeling of the merger history of dark matter haloes. 



LBG survey, we find that our estimated number densities are con- 
sistent with previous studies of LBGs in this redshift range. Over- 
all we obtain a mean redshift of 5 = 2.85 ± 0.34. From the 
data obtained we have shown evidence for the existence of galac- 
tic outflows with comparable offsets between emission and abs orp- 
tion lines as in previous studies (e.g. lPettini et al.ll200CI . 120021 and 
IShaplev et al. 20 03) 

We have further measured the clustering properties of the VLT 
VIMOS LBG sample. Based on the angular auto-correlation func- 
tion of the photometric LBG candidates, the real-space LBG cor- 
relation function, £,{r), is estimated to take the form of a dou- 
ble power-law, with a break at r^ ~ 1.5 ft^^Mpc. This is 
parametrised by a clustering length and slope below the break of 
ro,i = 3.19 ± 0.55 /i"^Mpc, 71 = 2.45 ± 0.15 and above the 
break of ro,2 = 4.37lg:5| /i"^Mpc, 72 = 1.61 ± 0.15. 

Assuming 7 = 1.8, the semi-projected LBG correlation 
function vjp{a) gives tq = 'i.(iJ^Q2A h~^yipc for the VLT 

= 4.2l°:J^ ft-^Mpc for the 



LBGs, slightly lower than ro 
Keck LBGs, and the combined VLT+Keck sample gives ro = 
3.98lo;i2 h'^Mpc. At r;, > 1 /i^^Mpc, the ^(r) estimates from 
w{6) and Wp{a) are therefore quite consistent. At rj, < 1 h~^Mpc 
the steeper power-law from the angular correlation function rises 
above the single power-law that best fits Wp, but the difference 
is only marginally statistically significant. These measurements of 
LBG clustering are broadly consistent with previo us measurements 
of the clus tering of LBGs at z ~ 3 made by lAdelberger et al.l 
(120031) and Ida Angela et al.l ( l2005ar) but lower than those made by 
iFoucaud et alj J2003h 

We then measured the redshift-space LBG auto-correlation 



function, £,{s). As expected, this presents a flatter slope at scales 
s < 8 h~^Mpc due to the effect of velocity errors, outflows and 
intrinsic velocity dispersions. Both the VLT and Keck samples re- 
quire total pairwise velocity dispersions in the range < w^ >i/2 = 
600 - 1000 kms"^ to fit ^(s) higher than the < w f >'^^'^ = 
400 kms~^ previously assumed dda Angela et al.l2005ah . The VLT 
and Keck samples' ^(s) results both imply an intrinsic pairwise 
velocity dispersion of ±400 kms~^ for a ^(r) model with ro — 
3.98 h'^Mpc and 7 = 1.8. A higher < w^ >^''^ will imply a 
higher infall parameter, P{z = 3), due to the degeneracy between 
these parameters. The high value of the velocity dispersion will 
also have an impact on our search for th e effects of star-forma tion 
feedback on the QSO Lyman-Q forest dCrighton et al.ll201C f) be- 
cause any sharp decrease in absorption near an LBG will tend to 
be smoothed away by this dispersion acting as an effective redshift 



We combine our LBG sample with that of lSteidel et al.l ( l2003h 
with the aim of measuring the infall parameter, I3{z = 3). Us- 
ing a single power-law with ro = 3.98 h~^Mpc and 7 — 1.8 as 
our model for the real space ^(r), our fits to our measurement of 
the LBG S,{a, n) from the combined data-set produce a best fitting 
infall parameter of /? = 0.48 ± 0.17. We find that this value is con- 
sistent with the P — 0.37 value expected in the standard ACDM 
cosmology. For this cosmology the value of the LBG bias implied 
from the galaxy dynamics is & = 2.06l|!|g3, again consistent with 
the value of 6 = 2.22 ± 0.16 measured from the amplitude of the 
LBG ^{r), assuming the standard cosmology. 

We have als o ma de the cosmological te st suggested by 
iHovle et alj ( I2OO2I) and Ida Angela et alj hoOSah and shown that 
the values of fij^ and I3{z = 3) derived from LBG redshift-space 
distortion are consistent with those derived by comparing the am- 
plitude of LBG clustering ai z — 3 from the combination of the 
measured 2dFGRS clustering amplitude and 13 at z = 0.1, using 
linear theory. Our measurement of /3{z — 3) is therefore consis- 
tent with what is expected from the gravitational growth rate pre- 
dicted by Einstein gravity in the standard cosmological model (see 
.Guzzo et al...2008.) . 

Finally, we have used the clustering amplitude measured for 
the LBGs to test simple models of clustering evolution. In particu- 
lar, we find that if the LBGs are long-lived then they could be the 
progenitors of low redshift L* spirals or early-type galaxies by the 
present day. 

The VLT LBG Survey is an ongoing project and we hope to 
double the survey area and LBG numbers by completion of the 
project. In combination with this work we are performing a survey 
of 2 ~ 3 QSOs in our LBG survey fields using the AAOmega 
instrument at the AAT. Bringing these two data-sets together will 
present a significant data resource for the study of the relationship 
between galaxies and the IGM at z ~ 3. 
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